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ABSTRACT 

We study the escape of Lya photons from Lya emitting galaxies (LAEs) and 
the overall galaxy population using a sample of 99 LAEs at 1.9 < z < 3.8 detected 
through integral-field spectroscopy of blank fields by the HETDEX Pilot Survey. 
For 89 LAEs with broad-band counterparts we measure UV luminosities and 
UV slopes, and estimate E(B — V) under the assumption of a constant intrinsic 
UV slope for LAEs. These quantities are used to estimate dust-corrected star 
formation rates (SFR). Comparison between the observed Lya luminosity and 
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that predicted by the dust-corrected SFR yields the Lya escape fraction. We also 
measure the Lya luminosity function and luminosity density {piya) at 2 < z < 4. 
Using this and other measurements from the literature at 0.3 < z < 7.7 we trace 
the redshift evolution of pL ya - We compare it to the expectations from the star- 
formation history of the universe and characterize the evolution of the Lya escape 
fraction of galaxies. LAEs at 2 < z < 4 selected down to a luminosity limit of 
L(Lya) > 3 - 6 x 10 42 erg s" 1 (0.25 to 0.5 L*), have a mean (E(B - V)) = 
0.13 ± 0.01, implying an attenuation of ~ 70% in the UV. They show a median 
UV uncorrected SFR =11 M yr~ 1 , dust-corrected SFR = 34 M yr _1 , and Lya 
equivalent widths (EWs) which are consistent with normal stellar populations. 
We measure a median Lya escape fraction of 29%, with a large scatter and values 
ranging from a few percent to 100%. The Lya escape fraction in LAEs correlates 
with E(B — V) in a way that is expected if Lya photons suffer from similar 
amounts of dust extinction as UV continuum photons. This result implies that a 
strong enhancement of the Lya EW with dust, due to a clumpy multi-phase ISM, 
is not a common process in LAEs at these redshifts. It also suggests that while 
in other galaxies Lya can be preferentially quenched by dust due to its scattering 
nature, this is not the case in LAEs. We find no evolution in the average dust 
content and Lya escape fraction of LAEs from z ~ 4 to 2. We see hints of a 
drop in the number density of LAEs from z ~ 4 to 2 in the redshift distribution 
and the Lya luminosity function, although larger samples are required to confirm 
this. The mean Lya escape fraction of the overall galaxy population decreases 
significantly from z ~ 6 to z ~ 2. Our results point towards a scenario in 
which star-forming galaxies build up significant amounts of dust in their ISM 
between z ~ 6 and 2, reducing their Lya escape fraction, with LAE selection 
preferentially detecting galaxies which have the highest escape fractions given 
their dust content. The fact that a large escape of Lya photons is reached by 
z ~ 6 implies that better constraints on this quantity at higher redshifts might 
detect re- ionization in a way that is uncoupled from the effects of dust. 

Subject headings: galaxies: evolution, ISM, luminosity function, dust, extinction 



1. Introduction 



Lya photons are produced in large amounts in star forming regions, therefore it was 
predicted nearly half a century ago that the Ly a emission line at 12 16 A s hould be a sign- 
post for star-forming galaxies at high redshift ([Partridge fc Peebles! 119671 ) . Actual obser- 
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vations of Lya emitti ng (LAE) galax ies at high redshift had to wait for the advent of 8- 
10m class telescopes (IHu et al.lll998l ). A little more than a decade has passed since their 
discovery, and thanks to a series of systematic surveys at optical and near-infrared wave- 
lengths, large samples of LAEs, usually containing from tens to a few hun dred objects, 



have 



1998 



jeen compiled over a wide range o 



Rhoads et al. 



2000; 



Kudritzki et al. 



6; 



redshifts from z ~ 2 to z ~ 7 (eg. 



2000 



iMalhotra fc Rhoads 



7; 



2002: 



Ouchi et al. 



Cowie fc Hu 
20031 



Gawiser et al 



2006a 



Aiiki et al. 



200 



Gronwall et al. 



200 



Ouchi et al. 



2008 



Nilsson et al. 



2009a 



Finkelstein et al. 



2009; 



Guaita et al.ll2010tlHayes et al. 



2010a 



Ono et al. 



2010 



Adams et al 



201ll ). Space based ultra-violet (UV) observations have also been us ed to study Lya emitting 



galaxies at lower redshifts, all the way down to the local universe (IKunth et al.lll998l . 12003 
Haves et alibood . boorl : ktek et al.llioosl : beharveng et al.ll2008l : ICowie et al.lboioh . 



The intrinsic production of both Lya and UV continuum photons in a galaxy is directly 
proportional to the number of ionizing photons produced by young sta rs, which is propor- 
tional to the star formation rate (SFR) ( ]Kennicuttlll998l ; ISchaererll2003l ). In practice, we do 
not expect the observed Lya luminosity of galaxies to correlate well with their SFR because 
the resonant nature of the n=l-0 transition in hydrogen makes the escape of Lya photons a 
non-trivial radiative process. 

In principle, the large number of scatterings suffered by a Lya photon before escaping 
the neutral medium of a galaxy increase its probability, with respect to that of contin- 
uum photons outside the resonance wavelength, of being absorbed by a dust grain. Hence, 
we would expect even small amounts of dust in a galaxy's inte r-stellar medium (ISM) to 



severely decrease th e equivalent width (EW) of the Lya line ( jHummer fc Kunaszl 11980 



Chariot fc Falllll993l ). In reality the situation is far more complicated, and it is not clear 
how the extinction suffered by Lya, and that suffer ed by continuum photons, relate. One 



scenario which has b een proposed by several authors (jNeufeldlll991t lHaiman fc Spaanslll999 



Hansen fc Ohlbood ) is the possible enhancement of the Lya EW due to the presence of 



very clumpy dust distribution in a multi-phase ISM. For this type of ISM geometry most of 
the dust lives in cold neutral clouds embedded in an ionized medium. In this scenario, Lya 
photons have a high probability of being scattered in the surfaces of these clouds, spending 
most of their time prior to escape in the inter-cloud medium and actually suffering less dust 
extinction than non-resonant radiation, which can penetrate into the clo uds where it has 



a high er chance of being absorbed or scattered by dust grains. Recently, iFinkelstein et al. 



( 120091 ) claimed that this process can simultaneously explain the Lya fluxes and continuum 
spectral energy distributions of many objects in their sample of LAEs at z ~ 4.5. 

At high redshift the Lya line can also be affected by scattering in the inter-galactic 
medium (IGM), as escaping Lya photons bluewards of the line center can be redshifted into 
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the resonance wavelength. This effect is particularly important at z > 5 as the density of 
neutral gas in the universe increases, but even at lower redshifts, when the universe is almost 
completely ionized, intervening Lya forest absorption can occur. To first order, the IGM 
transmission b lue-wards of Lya is ~ 90%, 70%, and 50% at z ~ 1.9, 3.0, and 3.8 respectively 
(jMadaulll995l ). In the naive case where the line profile escaping a galaxy is symmetric and 
centered at the Lya resonance, since only photons bluewards of the line are affected, we 
can expect attenuations of ~ 5%, 15%, and 25% on the emerging flux at these redshifts. In 
reality the process can be significantly different. While inflow of IGM gas onto ga laxies can 
introduce further attenuation red- wards of the line resonance (IDijkstra et al.ll2007l ). outflows 
in a galaxy's ISM can redshift the emerging spectrum so as to be completely unaffected 
by the IG M (IVerhamme et al.l 120081) . For example, in a sample of 11 LBGs and LAEs at 
z ~ 3 — 5, Verhamme et al. (boOST ) find no need to introduce IGM absorption to successfully 
fit the observed line profiles. This, combined with the inherent stochasticity of intervening 
absorption systems towards different lines of sight, makes Lya IGM attenuation corrections 
very difficult and uncertain. 



mgs 



The kinematics of the neutral gas inside a galaxy and in its imm ediate surround 
also play an important role regarding the escape of Lya photons (IVerhamme et al. 



2006 1; lUiiRstra et al.ll^uubl; Irlansen upll^uubl; lUiiRstra et al.l I2UUYI : ivernamme et al.l I2UU8I : 
Adams et alil2009l : lLaursen et alJboid : Izheng et alJboioh . Simply put, the velocity field of 



Diikstra et a 



2006 



Hansen O 



2006 



Diikstra et all 120071 : IVerhamme et all 120081 : 



the neutral gas has a strong influence on the emission line profile of the Lya line. Different 
combinations of geometry and velocity fields can "move" photons out of the resonance fre- 
quency either by blueshifting (typically due to in-falling gas) or redshifting (due to outflows) 
them, changing the number of scatterings photons experience before exiting the galaxy as 
well as their escape frequency. This process can affect the amount of dust extinction as well 
as the amount of potential IGM scattering those photons will suffer. 

No clear agreement is found in the literature regarding the amount of dust present in 
the ISM of Lya emitting galaxies. While most studies of narrow-band selected LAEs at 
z ~ 3 seem to indicate they are consistent with very low dust or dust-fre e stellar population s 
jGawiser et all l2006al . 120071 : iNilsson et all I2OO7I : bronwall et all I2OO7I : buchi et al.l £3), 
there have been recent results suggesting that the LAE population is more heterogeneous 



and includes more dusty and evolved galax ies, especially at lower redshifts (ILai et al.l 12008 
Nilsson et al1l2009al : iFinkelstein et al.ll2009h . 



We use a new sample of spectroscopically detected LAEs at 1.9 < z < 3.8 from the The 



Hobb y Eberly Telescope Dark Energy Experiment (HETDEX) Pilot Survey (lAdams et al. 



20 111 ) to investigate the shape of the UV continuum of LAEs, as well as the Lya luminosity 
function of these objects, and to address: 
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• The dust content of LAEs, parameterized by the dust reddening E(B — V), and its 
evolution with redshift. 

• The star- formation properties (SFR), the Lja escape fraction in LAEs, and its evo- 
lution with redshift. 

• The relation between the dust content and the escape fraction of Lya photons. 

• The relation between the dust extinction seen by continuum and resonant Lya photons. 

• The contribution of LAEs to the integrated star formation rate density at different 
redshifts. 

• The Lya escape fraction of the overall galaxy population and its evolution with redshift. 



These galaxies have been detected through wide integral-field spectroscopic mapping 
of blank fi elds, using the Visible Integral-field Replicable Unit Spectrograph Prototype 



f VIRUS-P. iHill et al 



2008). The Pilot Survey catalog of emission line galaxies is presented 



Adams et al. ^2011 ) . hereafter Paper I. The large redshift range spanned by our sample 



m 

allows us to check for any potential evolution in the above properties of LAEs. 

In §2 we describe the HETDEX Pilot Survey from which the sample of Lya emitting 
galaxies is drawn. In §3 we present our sample of LAEs along with their luminosities and 
redshift distribution. §4 presents our measurement of the UV continuum slope and deriva- 
tion of the amount of dust extinction present in these objects. Discussion of any potential 
evolution in the dust properties of LAEs is also in this section. We compare both uncorrected 
as well as dust-corrected SFRs derived from both UV and Lya in §5, where we also compute 
the escape fraction of Lyct photons and show how it depends on the amount of dust redden- 
ing. In §6 we present the Lya luminosity function and check for its possible evolution with 
redshift. We compare the integrated SFR density derived from the Lya luminosity function 
to that for the global galaxy population in §7. In this way we can assess the contribution of 
LAEs to the star-formation budget of the universe at these redshifts and estimate the Lyo; 
photon escape fraction for the overall galaxy population. Finally, we summarize our results 
and present our conclusions in §8. 



H„ 



Throughout the paper we adopt a standar d set of ACDM cosm ological parameters, 
= 70 km s^Mpc -1 , tt M = 0.3, and tt A = 0.7 dDunklev et aDbopgh . 
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2. The HETDEX Pilot Survey 

Ever since their discovery, the standard method for detecting and selecting LAEs has 
been through narrow-band imaging in a passband sampling the Lya line at a given redshift. 
The redshift range of these type of surveys is given by the width of the narrow-band filter 
used, and is typically of the order of Az = 0.1. Hence, these studies are limited to very narrow 
and specific redshift ranges. In terms of surveyed volumes this limitation is compensated by 
the large fields of view of currently available optical imagers which allow for large areas of 
the sky (~ ldeg 2 ) to be surveyed using this technique. 

An alternative technique, which has been attempted for detecting LAEs over the last 



very low resolution slit-less spectroscopy 


(Kurk et a 


. 2004 


Deharvene et al. 


200£ 


0, blind 


slit spectroscopy ( 


M 


artin & Sawicki 2004: ' 


Iran et al. 


20041 I 


lauch et al. 


2008; 


Sawicki et al. 


2008; Cassata et al. 


2011 


). or integral- field spectroscopv (Ivan Breukelen et al. 


2005 


• 



The success of this type of surveys has been variable. While early attempts to detect 
LAEs at z ~ 6 using slit spectroscopy failed t o do so, and cou ld only set upper limits to 



their number density ( jMartin fc Sawickil 12004 



Tran et al. 



2004), mo re recent attempts a t 



lower redshifts (2 < z < 6) like the ones by iRauch et al. and ICassata et al.l ( 120111) 



have produced large samples of objects. Similarly, an early attempt by iKurk et al.l (120041 ) 
to find LAEs at z — 6.5 using slit-less spectroscopy only yielded one detection, while more 
recently space-based UV slit-less spectroscopy with the GALEX telescope has allowed for 
the construction of a large sample of LAEs at z ~ 0.3 (IDeharveng et al.l 120081 ) . The only 
att empt to detect LAE s usin g integral-field spectroscopy previous to this work was done 
by Ivan Breukelen et al.l ( 120051 ). who used the Visible Multi-Object Spectrograph (VIMOS) 
integral field unit (IFU) on the Very Large Telescope (VLT) to build a sample of 18 LAEs at 
2.3 < z < 4.6 over an area of 1.44 arcmin 2 corresponding to the VIMOS IFU field-of-view. 

Although when doing spectroscopic searches for LAEs the wavelength range, and hence 
the redshift over which Lya can be detected, is tens of times larger than for narrow-band 
imaging, surveyed volumes have been typically small due to the small areas sampled by the 
slits on the sk y, or the small fields- o f-view of most integral field units. For example, the 
IFU survey by Ivan Breukelen et al.l ( 120051 ) only covered ~ 10 4 Mpc 3 becau s e of t he small 
area surveyed, while the z ~ 3.1 narrow-band survey by iGronwall et al.l (120071 ) covered 
~ 10 5 Mpc 3 over a very narrow range of Az = 0.04 because of the large 36' x 36' area which 
can be imaged with the MOSAIC-II camera. It is clear that the most efficient way of building 
large samples of LAEs would be to conduct spectroscopic searches over large areas of the 
sky. 
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HETDEX jHill et al.l l2008h will survey ~ 60 deg 2 of sk jj| using the Visible Integral- 
field Replicable Unit Spectrograph (VIRUS. Eill et al.lboioh . a wide field of view (16' x 16') 
integral field spectrograph currently being built for the 9.2m Hobby Eberly Telescope (HET). 
HETDEX will produce a sample of ~ 8 x 10 5 LAEs at 1.9 < z < 3.5 over a volume of 8.7 
Gpc 3 . The power-spectrum of the spatial distribution of these objects will be used to set 
a percent level constrain on the dark energy equation of state parameter w at these high 
redshifts ( JHill et al.ll2008l ). A prototype of the instrument, VIRUS-P, is currently the largest 
field-of-view IFU in existence, and has been used over the last 3 years to conduct a Pilot 
Survey for LAEs from which the sample used in this work is taken from (Paper I). The 
Pilot Survey, described below, samples the 1.9 < z < 3.8 range, and covers a volume of 
~ 10 6 Mp c 3 over an area of 16 9 arc min 2 . This v o lume is ten times larger than the one 



Ouchi et al 



covered in iGronwall et al. (120071) and iGuaita et al.l ( 120101 ). three times larger than the one 
covered by iNilsson et al.l ( )2009af ). and of comparable size to the one sampled at z = 3.1 by 



( 20081 ) but over an area 20 times smaller, exemplifying the power of integral 



field spectroscopy to search for emission line galaxies over large volumes. 

The HETDEX Pilot Survey obtained integral field spectroscopy over ~169.23 arcmin 2 
of blank sky in four extra-galactic fields (COSMOS: 71.6 arcmin 2 , GOODS-N: 35.5 arcmin 2 , 



2003; 


Drorv et al. 


2001; 


Pierre et al. 


2004) 



scope at McDonald Observatory. The goal of the survey is to conduct an unbiased search 
for spectroscopically-detected emission line galaxies over a wide range of redshifts. Although 
a powerful dataset itself, the Pilot Survey also provides a proof of concept and a crucial 
test-bench for the planned HETDEX survey. 

The observations and data reduction, as well as the detection and classification of emis- 
sion line galaxies, are presented in Paper I, and we refer the reader to it for a more detailed 
description of the survey design. Briefly, each field is mapped by a mosaic of 1.7' x 1.7' 
VIRUS-P pointings (27, 13, 16, and 4 pointings in COSMOS, GOODS-N, MUNICS-S2, 
and XMM-LSS respectively). The VIRUS-P IFU consists of an square array of 246 fibers, 
each 4.235" in diameter, sampling the field with a 1/3 filling factor. While a set of three 
dithered exposures covers the field-of-view almost completely, we observed each pointing at 
six dithered positions, ensuring complete coverage and improving the spatial sampling of 
the field and the astrometric accuracy of our detections. For each pointing, we obtained 
spectra at 1,476 (6 x 246) positions, with any point on the sky being typically sampled by 2 
overlapping fibers. Overall, the Pilot Survey consists of ~88,000 individual spectra over 169 



1 The actual HETDEX footprint corresponds to a 420 deg 2 area, but only 1/7 of the field will be covered 
by fibers 
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arcmin of blank sky. Each spectrum covers the 3600A-5800A wavelength range with ~ 5A 
FWHM resolution (a inst ~ 130 km s" 1 at 5000A). 

After the data are reduced and a ID flux-calibrated spectrum is extracted for each fiber 
position, we search the "blank" spectra for emission lines using an automated procedure 
(Paper I). Line detections are associated, when possible, with counterparts in broad-band 
images available for all four fields. The VIRUS-P wavelength range allows the detection of 
common strong emission lines present in star-forming galaxies such as Lya at 1.9 < z < 3.8, 
[OII]A3727 at z < 0.56, H/3 at z < 0.19, [OIII]A4959 at z < 0.17, [OIII]A5007 at z < 0.16, as 
well as typical AGN lines like CIVA1549 at 1.3 < z < 2.7, CIII]A1909 at 0.9 < z < 2.0, and 
MgIIA2798 at 0.3 < z < 1.1. 

Source classification is based on the presence of multiple spectral lines when available. In 
the case of single line detections, the spectral classification is considerably more challenging. 
For LAEs, only the Lya line appears in our wavelength range, so we expect single line 
detections for our objects of interest. Nevertheless, [Oil] emitters at 0.19 < z < 0.56 will 
also appear as single line detections in the VIRUS-P spectra. Even [Oil] emitters at z < 0.19 
that have unfavorable emission line ratios can appear as single line detections if H/3 and the 
[OIII] doublet are below the noise level. Our 5A FWHM spectral resolution is not high 
enough to resolve the [OII]A3727 doublet, so we cannot rely on the line profile to classify 
these objects. While galaxies detected in redder lines such as H/3 and [OIII]A5007 can also 
appear as single line detections depending on their redshifts and line ratios, the volume over 
which we sample these galaxies is ~ 400 times smaller than the volume over which we sample 
LAEs, and ~ 20 times smaller than the volume over which we sample [Oil] emitters. Hence, 
contamination from H/3 and [OIII] emitters is negligible. 

The classification of single line detections is thoroughly discussed in Paper I, and is 
based on an EW criterion, where objects showing rest-frame EW(Lya) > 20A are classified 
as LAEs (for 4 objects the EW > 20A criterion was bypassed due to the existence of further 
evidence pointing towards their LAE nature; see Paper I). This EW constraint effectively 
reduces the contamination from low-z interlopers to a negligible level. A total of 105 Lya 
detections are present in the Pilot Survey catalog presented in Paper I. Of these, 6 show 
X-ray counterparts indicating an AGN nature, leaving a final sample of 99 "normal" star- 
forming LAEs. In Paper I we also present a thorough assessment of the completeness and 
spurious source contamination in our catalog, based on simulated data. The completeness 
is used in §6 to estimate the Lya luminosity function. In our sample of LAEs we expect a 
4-10% contamination from spurious sources. The sample used in this work is presented in 
Table HJ 
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3. LAE Sample 

The 99 LAEs in the sample span a range in luminosities of log(LL ya ) = 42.42 — 44.03, 
and have a median luminosity of log(Li J/a ) = 43.03. Figure [1] shows the survey 5a limiting 
Lya luminosity as a function of redshift, together with the luminosities and redshifts of all 
LAEs in the sample. The depth of the observations is variable across the survey area and 
dependent on the observing conditions, the airmass at which the observations were taken, 
the Galactic dust extinction towards different fields, and the instrumental configuration. 
Colors in Figured] correspond to the fraction of the total surveyed area for which the spectra 
reaches the corresponding limit in luminosity. While VIRUS-P has its lower throughput 
in the blue end of the wavelength range, the smaller luminosity distance at lower redshifts 
compensates for this fact, providing a relatively flat luminosity limit throughout the entire 
redshift range. As mentioned above, detailed simulations quantifying the completeness and 
spurious detection ratio for the whole survey are presented in Paper I. A good understanding 
of the completeness of the survey is essential in order to calculate the Lya luminosity function. 
As shown in Paper I, the completeness at the 5a flux limit shown in Figure [T]is 33%, reaching 
50% at 5.5a and 90% at 7.5a. 

The redshift distribution of LAEs in our sample is shown in Figure [2] (errorbars show 
Poisson statistical uncertainties). The detected galaxies span a range in redshift of 2.079 < 
z < 3.745, with a median redshift of z = 2.811, properly sampling the 1.9 < z < 3.8 range 
over which they could be detected. Figure |2] also shows t he predicted re d shift distribution 



of LAEs in the Pilot Survey calculated by integrating the iGronwall et al.l ( 120071 ) luminosity 
function of narrow-band selected LAEs at z — 3.1 above the Pilot Survey flux limits shown 
in Figure [TJ and correcting for the survey completeness. The agreement is excellent at high 
redhsift (z > 3), but we observe a drop in the number of LAEs at lower redhsifts from 
what is predicted by a non-evolving lumin osity function. Re cent narrow-ba nd studies of 



z ~ 2 LAEs show hints for both an increase (IGuaita et al.ll2010l ) and decrease ( jNilsson et al. 



2009af ) of the LAE number density from z = 3 to z = 2. As stated by the authors themselves, 
neither of these studies probe a large enough volume to allow for a significant detection of 
the evolution in the LAE number density. In our surveyed volume, which is a few times 
larger than the volumes surveyed in those studies, we find some evidence for a decrease in 
the number density of LAEs from z ~ 4 to z ~ 2, although as discussed in §6, the statistical 
uncertainties remain too large to make a definitive statement. In any case, this type of 
evolution is expected if the escape fraction of Lya photons from galaxies decreases towards 
lower redshifts. In §7 we find evidence that this effect indeed occurs, which supports the 
observed drop in the LAE number density. 
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4. The UV slope of Lyman alpha emitters 



The UV continuum slope has been shown to be a powerful t ool for estimating the amo unt 



of dust extinction in star - forming galaxies i n the local universe (Meurer et al 



well as at high redshift ( iDaddi et al.l 12004 ; iBouwens et al.l 120091 ; iReddy et al 



995 



19991 ) as 



20101 ) . Direct 



observations of the ultra-violet spectral energy distribution of local star- forming galaxies have 
demonstrated that in the 1000A-3000A range, the y are very well described by a power-law 
spectrum of the form fx oc X 13 (ICalzetti et al.lll994l ). Differential dust extinction (reddening) 
makes the power-law slope correlate well with the amount of dust extinction in galaxies. 

Measuring the spectral slope of the UV SED of LAEs at 1.9 < z < 3.8 provides a direct 
measurement of their dust content, and its evolution with time. Knowledge of the amount of 
dust extinction in LAEs allows us to correct UV measured SFRs. An unbiased measurement 
of the SFR in these objects is not only important regarding the star-formation properties 
of these galaxies, but can also be used, together with the observed Lya luminosities, to 
estimate the escape fraction of Lyct photons from the ISM of these high redshift systems, 
and to study a possible evolution in this quantity. 



4.1. Measurement of the UV continuum Slopes, UV Luminosities, and Lya 

EWs 



We identify continuum counterparts of spectroscopically detected LAEs in our sample 
using publicly available broad-band optical images sampling the rest-frame UV SED of the 
objects. Multi-band aperture photometry is then used to measure their UV continuum slope 
(3) and UV luminosity as described below. 

For the purpose of measuring 8 we use the B, r + , i + , a nd z + images of the COSMOS 
and GOODS-N fields presented in ICapak et all J2004h and ICapak et all J2007h . the g, r, 
i, and z images of the XMM-LSS fie ld from the Canada-France-Hawaii Telescope Legacy 
Survey fCFHTLS. iMellier et al.llioosh Wl field, and the g, i, and z MUNICS-Deep images 
presented in Paper I. The identification and association with broad-band counterparts of 
our emission line detected objects makes use of a maximum likelihood algorithm which is 
described in detail in Paper I. Briefly, our astrometric uncertainty and the typical surface 
density of galaxies as a function of continuum brightness is used to identify the most likely 
broad-band counterpart for each LAE. The possibility of the emission-line source having no 
counterpart in the broad-band imaging is also considered. This can happen if the source is 
fainter than the sensitivity of the images or if the source is spurious. The no counterpart 
option is adopted if the probability exceeds that of all other possible counterparts. Only 9 
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out of 99 (9%) objects show no broad-band counterparts. This number is in good agreement 
with the 4-10% contamination expected from spurious detections in our LAE sample (Paper 
I, §3), although these objects could in principle be real and have very high EWs. In Paper 
I we showed that only one of them has significantly high EW given the limits that can 
be put using the depth of the broad-band images, while the large majority (8/9) show low 
signal-to-noise (S/N) detections (< 6.5) where the false detection ratio is the highest. For 
simplicity we omit these "no counterpart" sources from our analysis as we expect the large 
majority of them to be false detections. We also reject one other object (HPS-89) from our 
analysis because its broad-band counterpart photometry is catastrophically affected by a 
bright neighbor. 

Fluxes are measured in optimal 1.4xFWHM diameter color apertures, and scaled to 
total fluxes for each object using the ratio between V-band (g-band for MUNICS) fluxes 
measured in the color a perture and aperture-co rrected fluxes measured in a SExtracto r 
( jBertin fc Arnoutsl Il996l ) defined Kron aperture (IGawiser et al.l l2006al ; iBlanc et al.l 120081 ) . 
Any contribution from the measured Ly-alpha line to the broad-band fluxes is removed. 
While we do take into account IGM absorption when fitting for the UV slope, we decide to 
omit the U-band in the fits because in our redshift range the band includes the Lyman 912A 
break. Since IGM absorption is expected to be stochastic, an average line-of-sight correction 
might not apply to single objects. This leaves us with a B-band through z-band SED for 
each object. 

The approximate rest-frame wavelength range sampled by the above bands shifts from 
1500A-3000A at z = 1.9 to 900A-1900A at z = 3.8, so only the B-band is affected by Lyman 
forest absorption at the higher redshif t end of our range. Follow ing a similar methodology as 
that described in lMeurer et al.l ( 119991 ) and iReddy et al.l (120101 ) . we compute the UV contin- 
uum slope for each object by fitting the rest-frame UV SED with a power-law spectrum of 
the for m fx cx X 13 corr ected for IGM absorption at the corresponding redshift of each object 



using a iMadaul ( 119951 ) prescription. All available bands redwards the Lyman break are used 
to perform the fit. When an object is not detected in a particular band, we properly include 
the upper limit in flux given by the photometric uncertainty in the \ 2 minimization in order 
to not censor our data. The error in /3 is estimated from Monte Carlo simulations of 100 
realizations of the UV SED, where the fluxes in each band are varied within their photo- 
metric errors. This fitting also provides the UV luminosities at 1216A and 1500A which are 
used to estimate the Lya EW and the SFR, respectively. All these quantities are reported 
in Table HJ Figure |3] shows the UV continuum slope (3 as a function of redshift for the 89 
objects having continuum counterparts. 



In principle, the UV slope can depend not only on the amount of dust extinction, 
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but also on the age, metallicity, and initial mass function (IMF) giving rise to the stellar 
population. Extensive work can be found i n the literature regard i ng the se effects on the 
observed UV slope of star-forming galaxies. iLeitherer fc Heckmanl ( 119951 ) showed that for 
both instantaneous bursts and constant star formation synthetic stellar populations, changes 
of the order of A/3 = ±0.2 around a typical value of /3 ~ —2.3 are introduced by variations 
in age (1 Myr to 1 Gyr) and metallicity (0.1 Z Q to 2 Z Q ). They also find the UV slope 
to be la rgely insensitiv e to th e assumed IMF. This result is in good agreement with the 
work of iBouwens et al.l (120091 ). who demonstrate that th e UV slope dependence on dust is 
dominant over that on age, metallicty and IMF. They use lBruzual &: Charlotl ( 120031 ) models 
to show that changes by a factor of two in age and metallicity introduce changes of A/3 < 0.1. 



Schaerer &: Pellol (120051 ) also present a similar result. In Figure 1 of their paper it can be 



seen that for a range in ages of 1 Myr to 1 Gyr (encompassing the expected age range for 
LAEs), and metallicities between 1/50 Z Q and solar, both constant SFR population synthesis 
models, and single bursts younger than 10 Myr (time over which they can produce significant 
Lya emission) show variations in their UV slopes of A/3 = ±0.2 (AE(B — V) = 0.04). These 
systematics are smaller than the typical uncertainty in the measurement of (3 for LAEs in 
our sample. Therefore, by assuming a constant value for the intrinsic (dust-free) UV slope 
across our LAE sample, we can robustly estimate the amount of dust reddening directly 
from the observed values of /3 given an attenuation law. 

The right axis of Figure [3] shows the corresponding value of the reddening EiB — V), 
calculated assuming an int rinsic UV slope /3 n = —2.23 for a dust-free stellar population 



(IMeurer et al.lll999l ). and a lCalzetti et al.l (120001 ) extinction law. The value of /3 is derived 
from a fi t to the relation be tween the IR to UV ratio and /3 in a sample of lo cal starburst 
galax ies (IMeurer et al.lll999l). and r eproduces the observed relation at z ~ 2 (IReddy et al. 



20101 ) . Alt hough iReddv et all (I2010f) found young (< 100 Myr) z ~ 2 galaxies to lie slightly 
below the Meurer et al. (jl999 ) relation, and closer to that of Pettini et al. ( 1998 ). these 



two relations converge at low extinction and imply basically indistinguishable values for 
/3o. In order to take into account age and metallicity induced uncertainties in our error 
budget for the dust reddening, we sum in quadrature a systematic error of A/3 = ±0.2 
(AE(B — V) = 0.04) to the uncertainty in /3, and propagate it into the error in E(B — V). 
Measured values for the dust reddening and it associated uncertainty are reported in Table 

m 
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4.2. Dust Properties of LAEs and comparison to Previous Measurements 



Our LAEs show a mean UV continuum slope (0) = — 1.5±0.1 (formal error on the mean) 
corresponding to a mean (E(B — V)) = 0.16±0.02 (median E(B — V) = 0.13). The measured 
slopes span a relatively broad range of — 3 < < +2, with the large majority (83/89, 93%) 
of the objects having < (E(B -V) < 0.45). All objects with < O (i.e. bluer than the 
assumed intrinsic dust-free slope) are consistent with = 0q (i.e. E(B — V) = 0) within la. 

These slopes and reddenings are in rough agreement with previ ous measurem e nts o f 



For z = 2. 1 LAEs. IGuaita et a" 



(120101 ) 

)) and a 



narrow-band detected LAE broad-band colors. 

find a typical (B — R) ~ 0.2 (0 = —1.5 using equation 3 in iNilsson et al.l ( 12009a 
relatively u niform distribution in the —0.5 < (B — R) < 1 (—3.3 < < 0.7) range. Similarly, 
at z — 2.3, INilsson et al.l (j2009al ) find a median (B — V) = 0.14 c orresponding to = —1.4, 



Hayes et al.l ( ]2010ah used 



0-0.4 . At h igher 



with the bulk of their LAEs having -3.0 < < 2.0. At z = 2.2 
SED fitting to find that LAEs in their sample have a range in E(B — V) 
redshifts, usually lower levels of extinction are measured. At z — 3.l lNilsson et al.l ( 120071 ) find 
A v = 0.26t|j;# from fitting the stacked SED of 23 LAEs in the GOO DS-S field, correspondin g 
to E{B - V) = 0.06+ 

aol (assuming a Calzetti attenuation law). IVerhamme et al.l ( 120081 ). 
using Monte Carlo Lya rad iative transfer fittin g of the line profiles ofll,2~3 — 5 LBGs 



(8 of them also LAEs) fr omlTapken et al.l (120071) . find that the color excess spans a range of 
E(B - V) = 0.05 - 0.2. iGawiser et all J2006al ) report that the best-fit SED to the stacked 
optical photometry of z — 3.1 LAEs in their sample has Ay = OIq'o) corresponding to 
E(B - V) < 0.03. 



For compa rison, similar ranges in and E(B — V) as those seen her e have been measured 



for LBGs (e.g. iShaplev et al.l 1200 ll ; Erb et al.l 120061 : iReddv et al.ll2008l ). Figure 1 shows the 
E(B — V) distribution of LAEs in our sample, compared with that of UV c ontinuum s e lected 



gala xies at 1.9 < z < 2 7 (BX galaxies) and at 2.7 < z < 3.4 (LBGs) from Erb et al.l (120061 ) 
and lReddy et al.l ( 120081 ). The E(B — V) distributions for the continuum-selected galaxies are 
different from those presented in the original papers in that we have set all their E(B — V) < 
values to zero for proper comparison with our sample. It can be seen that both the shape 
and median value of the E(B — V) distribution of LAEs and BX/LBGs are relatively similar 
(medians are indicated by the dashed lines in Figure H|). This result, together with the 
fact that LAEs and BX/LBGs seem to ove rlap in the two-color BX/LBG selection diagram 
( IGuaita et al.l l2010l : IGawiser et al.l l2006al ). implies that both populations have relatively 
similar spectral continuum properties in the UV. Nevertheless, Figure H] shows an LAE 
distribution that is peaked at lower E(B — V) than the BX/LBG distributions, and also 
that Lya selection might allow for the inclusion of some highly reddened objects, although 
the reality of these red LAEs will be questioned in the next section. These galaxies, if real, 
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are excluded of UV-selected sample s by construction, s ince the color cuts in those selections 
reject object with E{B -V)> 0.5 JDaddi et al.ll2004l : blanc et al.ll2008h . 



The observed UV slopes imply that LAEs present low levels of dust extinction. One 
third (30/89) of the LAEs in our sample are consistent with being dust-free (E(B — V) = 0) 
to la, with the fraction going up to 60% within the 2a uncertainty. Still, a significant 
fraction of LAEs show non negligible amounts of dust. As will be shown in §5, dust in LAEs 
should not be neglected; doing so would strongly underestimate the SFR in these objects. 
Dust also plays a dominant role in setting the escape fraction of Lya photons, as we will 
discuss in §5.4. 



4.3. Evolution of the Dust Properties of LAEs 

At first sight, Figure [3] shows different behaviors in the dust-content distribution of 
LAEs at the high and low redshift ends of our sample. At z < 3 we see the emergence of 
a small population of LAEs (6/89) with red UV slopes (/3 > 0). This objects, if real, could 
represent an interesting population of dusty star forming galaxies in which some physical 
mechanisms allows for the escape of Lya photons. We have reasons to question the reality of 
these objects (see below). Furthermore, in this section we show that their presence does not 
affect the average properties of the LAE population which is dominated by UV-blue LAEs. 

To test for any evolution on the dust-content of LAEs with redshift we divided the 
sample in two redshift bins: low- z (1.9 < z < 2.8) and high- z (2.8 < z < 3.8). The division 
corresponds to the median redshift of the whole sample, and divides the survey volume in 
two roughly equal sub-volumes. The corresponding age of the universe at z = 3.8, 2.8, and 
1.9 is ~1.6, 2.3, and 3.4 Gyr. The median Lyct luminosity of the two sub-samples equals 
that of the whole sample (logLiya =43.0). 

The mean UV slopes of the high and low redshift samples are shown in Figure [3J Error- 
bars show the formal error on the mean and the standard deviation for each sample including 
~68% of the objects. The mean value of E(B — V) stays constant with values 0.16 ± 0.03 
and 0.17 ± 0.02 for the low and high redshift bins respectively. The scatter around these 
values is large and the means are statistically consistent with each other, and with the mean 
of the full sample. Therefore, we do not detect any significant evolution in the average UV 
slope and dust reddening of LAEs from z ~ 4 to z ~ 2. 

The lack of evolution in the mean dust-content of LAEs implies that this rare population 
of very high E(B — V) objects emerging at z < 3, if real, does not affect the average properties 
of the overall population due to their reduced number. The dust content of the bulk of the 
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LAE population remains relatively constant across the 2 < z < 4 range. 

Doubt regarding the validity of the UV slope measurements for these objects, and their 
classification as LAEs arises from looking at the distribution of rest-frame Lya EWs for 
our sample. Figure [5] shows the EW distribution for both UV-blue (/3 < 0) and UV-red 
((3 > 0) LAEs together with the one for the whole sample. The Lya EW is measured as 
described in §4.1, and hence can differ from the values presented in Paper I. It is evident from 
Figure that UV-blue LAEs dominate the overall population since they present a practically 
indistinguishable EW distribution (well fitted by an exponential with an e-folding parameter 
wo = 74 ± 7) from that of the full sample (wo = 77 ± 7). UV-red LAEs on the other hand, 
in addition to being rare in numbers, present a very different distribution in rest-frame EW, 
characterized by the presence of many extremely high EW (> 500A) objects. Two of these 
UV-red objects are in the MUNICS field where we lack deep X-ray data to reject AGNs from 
our sample (one of these sources shows significantly extended Lya emission and is a good 
candidate for an extended Lya nebula, or Lya Blob as discussed in Paper I). The remaining 
four objects have low association probabilities (< 0.6) with their broad-band counterparts, 
casting doubt on the validity of our UV slope and EW measurements for these objects. 
Further follow-up observations are necessary to confirm the nature of these detections. 

If real, these UV-red LAEs would have an extreme nature, being very dusty and highly 
star-forming. We remove these six objects from all the subsequent analysis, and for the rest 
of the paper we focus only on the results regarding the dominant UV-blue LAE population. 
It must be kept in mind that if these objects happen to be real LAEs, no strong evidence 
for a bimodality in the dust content or SER of LAEs is found in our data. The /3 > cut 
used to separate this population of objects is solely based on the fact that (3 > objects 
are absent at z > 3 in our sample. After removing these objects from our sample, we find 
a mean dust-reddening for LAEs of {E(B — V)) = 0.13 ± 0.01, corresponding to an average 
dust attenuation of ~ 70% at 1500A. 



5. UV versus Ly« SFRs and the Escape Fraction of Ly« Photons 

In this section we use the dust extinction values derived from the UV continuum slope 
in the previous section to estimate the dust-corrected SFR of LAEs in our sample. A 
comparison between the observed Lya luminosity and the intrinsic Lya luminosity implied 
by the dust-corrected SFR allows estimation of the escape fraction of Lyct photons from 
these galaxies. Throughout this analysis we have decided to neglect the effects of the IGM. 
As stated in §1, at these redshifts we expect attenuations for Lyct of no more than 5-25%, 
which is within our typical uncertainty for the Lya luminosity. Furthermore, if outflows are 
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common in LAEs, as many lines of evidence suggest, then IGM scattering at these redshifts 
may become even less important as most Lya photons leave galaxies red-shifted from the 
resonance wavelength (see discussion and references in §1). We start by comparing the 
observed (not corrected for dust) SFRs derived from the UV and Lya luminosities, then 
introduce the dust corrections, and finally estimate the escape fraction of Lya photons and 
study how it relates to the amount of dust reddening. 



5.1. Estimation of the Star Formation Rate and the observed SFR(Lya) to 

SFR(UV) ratio. 



The UV monochromatic luminosity at 1500A (-£^1500) for each object is taken fr o m the 
fits described in §4. In order to calculate the SFR we use a standard iKennicuttl (119981 ) 



conversion 



SFR(UV) [Moyr" 1 ] = 1.4 x KT^isooferg s^Hz" 



which assumes a Salpeter IMF wi th mass limits 0. 1 to 100 M Q . The Lya derived SFRs 
were calculated using the standard IKennicuttl (119981 ) conversion factor fo r Ha and assuming 



the intrinsic Lyq to Ha rat io of 8.7 from Case B recombination theory (IBrocklehurstl 11971 
Osterbrock fc Ferlandll2006l ). so 



SFR(Lya) [M^ 1 ] = 7.9 x l - 42 ^[erg s" 1 ] 



(2) 



Figure [6] shows SFR(Lya) versus SFR{UV) for our 83 objects. Without accounting for 
dust we measure median SFRs of 11 M Q yr _1 and 10 M yr _1 from UV continuum and Lya 
respectively. Although these agree with what is typically quoted for LAEs in the literature, 
we consider them to be underestimated by roughly a factor of ~ 3 — 4 because of the lack of 
a dust extinction correction. 

We observe a median ratio SFR(Lya)/ SFR{UV) = 0.83. Single objects present a large 
scatter around the median, with values ranging from 0.2 to 5.9. Since the UV SFR conversion 
factor is valid for galaxies with constant star formation over 100 Myr or more, while the 
one for Lya is valid at much younger ages of ~ 10 Myr (IKennicuttl Il998l ) . young galaxies 
can have intrinsic SFR(Lya) to SFR(UV) ratios higher than unity. The dashed lines in 
Figure [6] show the allowed range for dust-free constant star-formation st ellar popu lations 



with metallicities from 1/50 Z Q to solar and ages from 1 Myr to 1 Gyr from ISchaererl (120031 ). 



A Lya escape fraction of less than unity can push objects above this range. All the objects 
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in our sample show SFR(Lya) to SFR(UV) ratios (or roughly equivalently Lya EWs), 
which are consistent within la with those of normal stellar populations (i.e lower than ~ 4). 

The observed median ratio between these two quantiti es is in rough agree ment with 
previous measurements found in the literature. At z ~ 2.1, iGuaita et all (120101 ) measures 
a mean SFR(Lya) to SFR(UV) ratio of 0.66 for nar row-band sel ected LAEs, consistent 
with the 0.56 value measured by iNilsson et all ( l2009a[ ) at z ~ 2.3. iGronwall et all (120071 ) 
find LAEs at z ~ 3.1 to span a similar range in the SFR(Lya)-SFR(UV) plane as the 
one observed here, and while they quote a mean ratio o f 0.33, a revised value of ~ 1 is 
actually a better estimate for their samp le B Ouchi et al. ( 2008) m easures a ratio of 1.2 in 
their z ~ 3.1 sample of LAEs. Recently iDijkstra fc Westra I ( 120101 ) conducted a statistical 
study of the relation between these two quantities. Compiling a number of LAE samples at 
3-0 < z < 6.5, they find 68% of LAEs to show SFR(Lya) / SFR(UV) = 0.9+J|, in agreement 
with our observations. 

There is reason to expect evolution in SFR(Lya)/ SFR(UV) with redshift. First, if 
the dust content of galaxies changes with redshift, and Lya and UV photons suffer different 
amounts of extinction, we should see a redshift dependence in the ratio. Also, if the Lya 
line suffer from significant IGM absorption, the dependence of the IGM opacity with redshift 
should affect the SFR(Lya) to SFR{UV) ratio. In Figure [7] we present the SFR(Lya) to 
SFR(UV) ratio, as well as the rest-frame Lya EW as a function of redshift. While these 
two quantities are roughly equivalent, SFR(UV) is calculated from the UV monochromatic 
luminosity at 1500A, while the EW uses the monochromatic luminosity at 1216A, therefore 
the ratio between them has a mild dependence on the UV slope. Because of this dependence, 
we chose to present both quantities in Figures [7] and [HJ 

Over the 2 < z < 4 range, we do not observe evolution at a significant level in the Lya 
EW or the ratio between the Lya and UV SFRs. For our low and high redshift bins we 
measure median EWs of 87 ± 63A and 53 ± 26A respectively (median absolute deviation 
errors). A Kolgomorov Smirnof (KS) test to the cumulative EW distributions for the low 
and high redshift sub-samples allows the hypothesis of them being drawn from the same 
parent distribution to 2a. In terms of SFR(Lya)/ SFR(UV), the measured median ratios 
are 1.1 and 0.7 for the low and high redshift sub-samples. The fact that we do not observe a 
significant decrease in the typical EW of LAEs supports our assumption of neglecting IGM 
absorption in our analysis. 

We also analyze the relation between SFR(Lya)/ SFR{UV) and the dust reddening 
E(B — V) derived from the UV slope. If UV and Lya photons suffer from similar amounts 



2 Caryl Gronwall, private communication 
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of extinction, the above ratio should be independent of the amount of dust present in the 
galaxy. This is indeed the case for our LAEs, as can be seen in Figure where the relation 
for the two quantities (as well as that between EW and E(B — V)) is shown. Throughout 
the entire range < E(B — V) < 0.45 the ratio between Lya and UV derived SFRs stays 
flat with objects scattered around the median value. A similar behavior is seen for the EW. 



5.2. Dust Corrected SFRs and Estimation of the Lya Escape Fraction 

We no w correct the UV lum inosity of our objects using the values of E(B — V) estimated 



in §4 and a ICalzetti et all ( 120001 ) attenuation law. This approach provides a better estimate 
of the true SFR in the galaxies. Figure [9] shows a comparison between the dust-corrected 
SFR(UV) corr = SFR(UV) x io^iboo e(b-v)) ? and the uncorrec ted SFR(Lya). Error- 
bars include the uncertainty in the dust correction which has been propagated from the 
uncertainty in the measurement of the UV continuum slope /3. Note that the axes in Figure 
HI are different from those in Figure LAEs in our sample have a median dust-corrected 
SFR(UV) = 34 M Q yr _1 , a factor of ~ 3 higher than the uncorrected median value, and 
show intrinsic SFRs ranging from 1 to 1500 M yr" 



-l 



The escape fraction of Lya photons is given by the ratio between the Lya derived SFR 
and the extinction corrected UV SFR. 



n i T \ L{Lya) observed 

Jese{Lya) = — r 

l J \l-'y ) intrinsic 



SFR(Lya) 

SFR(UV) x io(o.4fcisoo£(B-v)) 

(3) 

Measured values of f esc (Lya) are reported in TableHJ Figure [TUI presents the Lya escape 
fraction of our LAEs as a function of redshift (black circles). A broad range in the escape 
fraction (2% to 100%) is observed. LAEs in our sample show a median escape fraction 
fesc(Lya) = 0.29 ± 0.04, and a mean escape fraction (f esc (Lya)) = 0.55 ± 0.08 (formal error 
on the mean). All objects showing f esc (Lya) > 1 are consistent with f esc (Lya) = 1 to 1.5cx. 



A recent study by iHayes et all ( l2010al ) used a pair of optical and NIR narrow-band 
filters to sample the Lya and Ha lines over the same volume. By comparing the Lya and 
Ha luminosities of a sample of 38 LAEs at z = 2.2, they derived a lower limit of 0.32 for 
the average Lya escape fraction of LAEs, a value consistent with our measured average. 
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Another estimation of the Lya escape fraction was done by IVerhamme et al.l (120081 ) using 
an independent method on their spectroscopic sample of 11 high-z galaxies (8 of which are 
LAEs). Fitting the Lya emission line velocity profiles using Monte Carlo radiative transfer 
simulations yielded best-fit values for f esc (Lya) between 0.02 and 1, with a median value 
of 0.17, in good agreement with our observed median value. The agreement between these 
three independent estimations using a different set of techniques is encouraging. 

A median (mean) escape fraction of ~ 20% (~ 50%) is one order of m agnitude higher 
than that adopted in the semi-analytical models of iLe Delliou et al.l ( 120051 ). in which a 2% 
escape fraction combined with a top-heavy IMF is used to match the output of the models 
to the observed Lya and UV luminosity function of LAEs at different redshifts. We observe 
a much larger escape fraction, and our measured EWs can be explained by standard stellar 
populations with normal IMFs. Also, the large scatter seen in Figure [TD] implies that using 
a single value of f esc (Lya) to model the LAE galaxy population is not a realistic approach. 

It is important to remark that estimating the escape fraction directly from the observed 
SFR(Lya)/ SFR{UV) ratio, by assuming LAEs are dust-free galaxies, would imply a sig- 
nificant overestimation of its val ue. For examp l e, the best-fit SED to the stacked optical 
photometry of z — 3.1 LAEs in iGawiser et al.l ( 2006al ). which has Ay = OIq'q, implies a 
best-fit escape fraction of 0.8 (although the uncertainty in the fit allows for a escape fractio n 
> 0.2, in agreement w ith o ur results). S i milarl y, the ratios measured by lOuchi et al.l (120081 ). 



Nilsson et al.l (j2009al ). and lGuaita et al.l (120101 ) imply escape fractions in the 0.5 to 1.0 range 



if dust is not considered. As discussed above, if we were to completely neglect dust extinction 
we, would measure a value of 0.87 for our sample. 



5.3. Evolution of the Lya Escape Fraction in LAEs 



No significant evolutionary trend is present across the 1.9 < z < 3.8 range of our objects 
in Figure [TUj where the median escape fractions for the two 1.9 < z < 2.8 and 2.8 < z < 3.8 
redshift bins (red open stars in Figure [TUj) are consistent with the median for the whole 
sample. In order to investigate if the Lya escape fraction of LAEs evolves over a larger 
baseline in cosmic time, we also show results found in the literature at a lower redshift. 
At higher redshifts the Lya escape fraction for L AEs remains poo rly constrained, although 
attempts to measure it exist in the literature (eg. lOno et al.ll2010l ) 



At lo w redshift lAtek et al.l (120091) performed optical spectroscopy on a sample of z ~ 
0.3 LAEs ( Deharveng et alT2008 ) and used the Ha luminosity, in combination with dust 
extinctions derived from the Balmer decrement, to estimate the Lya escape fraction of these 
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objects. A similar range in the escape fraction is observed for z ~ 0.3 and 2 < z < 4 LAEs, 
with the former showing values ranging from 0.03 to 1, implying that there has not been 
significant evolution in f Pfsr { Lya) over t he ~ 8 Gyr from z ~ 3 to z ~ 0.3. At very high 



redshifts (z = 5.7 and 6.6) lOno et al.l (120101 ) has estimated the Lya escape fraction of a 



sample of a few hundred narrow-band selected LAEs using a similar method to the one used 
here, except that their intrinsic SFRs were measured by SED fitting. Their escape fractions 
are consistent with our measured values at 2 < z < 4, although their error-bars are large. 
Therefore we detect no significant evolution in f esc (Lya) over the 0.3 < z < 6.6 range. 

This lack of evolution in the Lya escape fraction of LAEs must be interpreted with 
caution, since nothing ensures that the LAE selection technique recovers the same galaxy 
populations at these very distant epochs in the universe. Furthermore, since the selection 
is based on the strength of the Lya line relative to the the underlying continuum (i.e. the 
EW of the line), the technique will tend to favor galaxies with high Lya escape fractions, as 
long as they satisfy the brightness cut of the survey, at any redshift. Therefore, the lack of 
evolution in f esc (Lya) cannot be interpreted as constancy in the physical conditions in the 
ISM of these galaxies. For example, while at low redshift the escape fraction is most likely 
dominated by dust absorption, at z ~ 6 it is most likely dominated by IGM attenuation. 



5.4. The Relation between f esc (Lya) and Dust 



As discussed in §1, a major subject of debate regarding the escape of Lya photons from 
star-forming galaxies is the role played by dust. It is not clear whether the resonant nature 
of the transition produces Lya photons to be extincted more, less, or in the same amount as 
continuum photons outside the resonance. For example, while Lya photons should originate 
in the same regions as Ha photons, we have no reason to expect the extinction seen by Lya 
photons to follow the nebular extinction relation EiB — V) s t ar s = 0.4A E(B — V)ga.<i s een for 
non-resonant hydrogen transitions in star- forming galaxies at z = ( ICalzettil Il997l ). since 
resonant scatter makes the optical paths seen by Lya completely different from the one seen 
by other lines like Ha or H/3. Furthermore, it has not been established if the above relation 
holds at high redshift or not. 

In order to test this issue we parameterize the ratio between the dust opacity seen by 
Lya and that which cont inuum photons at the sa me wavelength would see in the absence of 
the transition. Following iFinkelstein et al.l (120081 ) . we adop t the parameter q = TL ya /rx = i2i6, 
where T\ = k\E(B — V)/1.086 and k\ is assumed to be a ICalzetti et al.l (120001 ) dust atten- 
uation law. EiB — V) is always taken to be the stellar color excess derived from the UV 
slope. 
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A value of q ~ implies that Lya photons suffer very little extinction by dust, as is 
expected in an extremely clumpy multi-phase ISM. Large values (q ^> 1) represent cases in 
which scattering of Lya photons introduces a strong increase in the dust attenuation as is 
expected in a more homogeneous ISM. As discussed in §1, not only the structure of the ISM 
determines the value of q, but also the kinematics of the ISM, since favorable configurations 
(eg. an expanding shell of neutral gas which allows backscattering) can reduce the amount 
of dust extinction seen by Lya photons. All these processes are coupled in determining the 
value of q, and discriminating between them requires a joint analysis of the UV, Lya, and 
Ha luminosities, the dust extinction either from the shape of the UV continuum or from 
Balmer decrements, and the profiles of the latter two emission lines. Until such data exist, 
interpretation is difficult, but we can still gain insights about the escape of Lya photons and 
the dust properties of LAEs from the measured value of q. 

In Figure [TT1 we present the Lya escape fraction versus the dust reddening E(B — V) 
for our sample of LAEs. A clear correlation between the escape fraction and the amount of 
dust extinction is seen. Also shown are the expected correlations for different values of the 
parameter q. The LAE population falls along the q = 1 relation. The median for the whole 
sample is q = 0.99 ± 0.44 (median absolute deviation error), implying that in most LAEs 
the Lya emission suffers a very similar amount of dust extinction to that experienced by UV 
continuum light. 



Our results show good agreement with those of lHayes et al.l (l2010al ). In their work, 5 
out of 6 LAEs at z — 2.2 detected in both Lya and Ha show escape fractions and dust 
reddenings consistent with the q = 1 relation. The same is true for the large majority of 
their LAEs with no Ha detections for which they could only provide lower limits to the 
escape fraction. Our LAEs at 2 < z < 4 also occ upy the same reg ion in the E(B — V) vs. 



f esc (Lya) plane as the low redshift LAEs (z ~ 0.3 lAtek et al.ll2009l ) shown as green triangles 
in Figure [TTJ This implies that not only the Lya escape fraction in LAEs does not evolve 
with redshift as shown in §5.3, but its dependence on the dust content of the ISM remains 
the same from z = 4 to 0.3. 

The LAE selection criteria imply that these objects are chosen to be the galaxies with 
the largest Lya escape fractions given their Lya luminosities and dust content at any redshift. 
Most likely, a combination of ISM geometry and kinematics favors the escape of Lya photons 
in these galaxies as compared to the common galaxy population at any redshift. Hence, 
when examining Figure [11] we should think of LA Es as the upp e r eny elope of the escape 



fraction distribution at any E(B-V). For example, iKornei et al.l (120101 ) found that LBGs 
with Lya emission typically lie below the q — 1 relation. In their work they parameterized 
the difference in the extinction seen by Lya and continuum photons by the "relative escape 
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fraction" , f eS c,reh which relates to q by the following relation 

■Xd ) t A \ 

9 = 0Ak x=1216 E(B - V) 1 J 

They find LBGs to have (f esC! rei) — 0.27 (which does not include LBGs showing Lya in 
absorption). We present this relation for LBGs as a red line in Figure [TTJ This finding 
supports our previous point, namely that LAEs are the upper envelope to the overall galaxy 
population in the E(B — V) vs. f esc (Lya) plane. 

Our result should not be interpreted as evidence against the existence of a clumpy 
multiphase ISM in LAEs, since in the presence of a completely homogeneous ISM we expect 
q > 1. Nevertheless, our result indicates that either a clumpy ISM, a favorable kinematic 
configuration of the ISM, or a combination of both, can reduce the amount of dust seen by 
scattering Lya photons only up to the point where they encounter the same level of dust 
opacity as the continuum. Since LAEs by definition will be the galaxies with the largest Lya 
escape fractions at any value of E(B — V), the absence of a significant number of points at 
low values of q in Figure [TTJ suggests that enhancement of the Lya EWs due to dumpiness 
in the ISM is not a common process in galaxies. 



6. The Lya Luminosity Function 



It has been well established in the literat ure that the Lya luminosity function does 
not show significa nt evolution from z = 3 to 6 (IShimasaku et al.l 120061 ; iTapken et al.l 12006 



Ouchi et al.l 120081 ). On the other hand, a strong decrease of roughly one order of magnitud e 
is seen in the abundance of LAEs from z ~ 3 down to z ~ 0.3 ( iDeharveng et al.l 120081 ) . 
At what point in cosmic time this decrease starts to take place, and how well it traces 
the SFR history of the universe, is unknown. Recently there ha ye been reports of po ssible 
evolution in the number density of LAEs between z ~ 3 and z ~ 2 ( INilsson et al.ll2009al ). but, 
as stated by the same authors, these result s might be aff e cted b y cosmic sample variance 
over the surveyed volumes. Furthermore, ICassata et al.l (120111 ) find no evolution in the 
luminosity function across these epochs in their sample of spectroscopically detected LAEs. 
The existence of evolution in the luminosity function (or number density) of LAEs down to 
these lower redshifts is still a subject of debate. 

By examining the redshift distribution of our LAEs (Figure [2]), we found indications that 
their number density might be decreasing towards lower redshifts (z < 3) in our sample (§3). 
In this section we measure the Lya luminosity function of LAEs, and study any potential 



-23 - 



evolution down to z ~ 2. We restrict the measurement of the luminosity function to LAEs 
in the COSMOS and GOODS-N fields, which account for 81% (80/99) of the total sample. 
Both MUNICS and XMM-LSS lack deep X-ray data comparable to the one available in 
COSMOS and GOODS-N, so it is not possible to identify AGNs in those fields. 



6.1. Measurement of the Luminosity Function 

To measure the lum inosity function we adopt a 1/V max formalism similar to the one used 



by ICassata et al.l (120111 ). We compute the volume density of objects in bins of Alog(L) = 
0.125 dex, as the sum of the inverse of the maximum volumes over which each object in 
the luminosity bin could have been detected in our survey. As discussed in §3, the depth of 
the observations is variable across the surveyed area. The whole survey covers 169 arcmin 2 , 
corresponding to 60 VIRUS-P pointings. Each pointing was covered by six dithered observa- 
tions, which accounts for 360 independent observations each reaching different depths. The 
noise spectrum for each IFU fiber in each of these observations is an output of our data 
processing pipeline, and can be directly translated into an effective line luminosity limit for 
Lya at each redshift (see Figure [TJ . 

For each object, V max is given by 



V m ax ^ V max i (5) 



where V max ^ is the integral of the co-moving volume over all the redshifts at which the 
object could have been detected (i.e. where the luminosity limit is lower than the objects 
luminosity) for each observation i. Summing over the inverse of V max for all objects in each 
luminosity bin then yields the luminosity function shown as open black circles in Figure [T2j 

As mentioned briefly in §3 and discussed extensively in Paper I, the effects of incom- 
pleteness are important over all luminosities in our survey. Completeness is a direct function 
of the S/N at which the emission line is detected in our spectra. Since the limiting luminosity 
is not constant at all redshifts (Figure [TJ, objects of the same luminosity can be detected with 
high significance, and hence high completeness, at certain redshifts and with low significance 
and low completeness at others. This is different than, for example, imaging surveys where 
objects are detected in a narrow redshift range, and the S/N is close to a unique function 
of the luminosity. In those cases, incompleteness becomes only important in low luminosity 
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bins, where the objects flux approach the depth of the images^. In our case, we must account 
for incompleteness over the whole luminosity range if we want to get a proper estimate of 
the luminosity function. 

In Paper I we present a detailed completeness analysis of our survey based on simulations 
of the recovery fraction of synthetic emission lines at different S/N in our spectra. Using these 
recovery fractions we correct the observed Lya luminosity function calculated as described 
above. The resulting completeness-corrected luminosity function is shown by the red filled 
circles in Figure [T2"l Error-bars shows Poisson uncertainties only, and co rrespond to a lower 
limit for the error since they do not include cosmic variance, although lOuchi et al.l (120081 ) 
show that for volumes such as the one surveyed here (~ 10 6 Mpc 3 ), cosmic sample variance 
uncertainties are not significantly higher than Poisson errors. 



We fit the observed Lya luminosity function using a ISchechterl (119761 ) function of the 
following form 



(p(L)dL = (p*(L/L*) a exp(-L/L*) d(L/L* 



(6) 



Since the depth of our observations (~5x 10 erg s cm in line flux) is somewhat 
limited, we do not consider our data to be sufficiently deep to constrain the faint-end slope 
(a) of the luminosity function. We consider the be st available constrain t on a to come 
from the spectroscopic survey recently performed by ICassata et al.l (120111 ) . They measure 
a ~ —1.7 using a survey which reaches more than one order of magnitude deeper than ours 
in terms of limiting line flux (~ 1.5 x 10~ 18 erg s -1 cm~ 2 ). We take their measured a as 
our fixed fiducial value for the faint-end slope of the luminosity function, but also report 
results assuming a = —1.5, since that i s the value most commonly used in the literature 
(IGronwall et al.l 120071 ; lOuchi et al.l 120081 ). Our results do not depend significantly on the 
assumed value of a. 



The best-fit Schechter luminosity function (a = —1.7) is shown by the solid red line in 
Figure IT2"l and 1, 2, and 3a confidence limits for the parameters are shown in Figure IT51 The 
best fit parameters for a = —1.7 and -1.5 are reported in the first two rows of Table [2J 



In reality, incompleteness in narrow-band emission line surveys is more complicated than this because 
of the non top-hat s hape of narrow-band filters, and shows a dependence with the redshift of the sources; 
see the discussion in IGronwall et al.l (|2007 ). 
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6.2. Comparison with Previous Measurements 
Figure [12] also shows the Lya luminosity functions measured by several authors at sim- 



nosity funct 


ions of 


van Breukelen et al. 


Haves et al. 


( 


2010a 


), and 


Cassata et al. 



(2005); iGronwall et all (120071 ): lOuchi et all (12008( 1: 



(120111 ) measured at 2.3 < z < 4.6, z = 3.1, z 



2.2, and 1.95 < z < 3 res pectively, agree with our observed values to within 



3.1, 
; la 



(Poisson) at all luminosities. The lHayes et al.l (j2010al ) measurement shows better agreement 



with our data at the bright end of the luminosity function. This is in fact surprising, as their 
measurement was performed over a smaller volume (5.4 x 10 3 Mpc 3 ) and a fainter range in 
luminosities (2 x 10 41 — 5 x 10 42 erg s _1 ) than the other mentioned works. 

Our best-fit Schechter function appears to be flatter than previous measurements over 
a similar range in luminosities. Figure [TBI shows that thi s difference i s becau se we derive a 
higher L* and a lower <p* than previous authors (except lHayes et al.l (l2010al ) who found a 
very similar value for 4>* but a larger value for L * ). The best-fit pa r ameters measured b y 



van Breukelen et al l (120051 ): IGronwall et al.l (120071 ) : lQuchi et al.l (120081 ) : lHaves et all (l2010ah . 
and Cassata et al. ( 20 111 ) fall within our 2a confidence contour. This last work is the only 
one of the mentioned luminosity function measurements in which a = —1.7. For all the 
ot her measurements, th e faint -end slope was either assumed or measured to be —1.5 except 
for Ivan Breukelen et al.l (120051 ) who used —1.6. For better comparison, Figure [TBI also shows 
uncertainty contours for our fit assuming a = —1.5 (dotted contours). As mentioned above, 
the value of a does not change our results in any significant way. 

The L* and (f>* parameters are strongly correlated with each other, so the 2a discrepancy 
with previous measurements is not surprising as it follows the sense of the correlation. Most 
importantly, we survey a very large volume and hence are able to find rare high luminosity 
objects. The luminosity functions derived in these studies typically stop at ~ 10 43 erg s _1 , 
while we find objects up to three times brighter luminosities. If we fit a Schechter function to 
only bins with L(Lya) < 10 43 erg s" 1 , we obtain the luminosity function shown as a dashed 
red line in Figure [T2J which is in much better agreement with previous measurements (black 
star in Figure [TBI and third row in Table [2]). 



6.3. Evolution of the Lyo Luminosity Function 

As mentioned above, evidence suggests that the Lya luminosity function does not signif- 
icantly evolve between z ~ 3 and z ~ 6. While at the high end of this redshift range [z >5) 
IGM absorption might become important and the lack of evolution might imply an increase 
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in the intrinsic Lya luminosity function (ICassata et al.ll201ll ). at least between z ~ 4 and 



z ~ 3 the lack of intrinsic evolut ion seems well supported as changes in IGM transmission 



are negligible (jOuchi et al.l 120081 ) . We can extend these studies to lower redshifts and ask: 



Does the luminosity function show any significant evolution down to z ~ 2? 

To test for possible evolution, we have again divided our sample in the two redshift 
bins defined in §4, one at low-redshift (1.9 < z < 2.8), and another one at high-redshift 
(2.8 < z < 3.8). We measure the luminosity function in each of these bins independently. 
The results are shown in Figure H2J best fit parameters are presented in Table [2J and la 
confidence limits are shown in Figure [TBI At L(Lya) < 10 43 erg s _1 , where cosmic variance 
is lower than at higher luminosities, the low-z luminosity function seems to be systematically 
lower than the high-z luminosity function by a factor about ~ 2, in agreement with what we 
observed in §3 when comparing the redshift distribution of our objects to the predictions for 
a non evolving luminosity function. Still, both the high-z and low- z luminosity functions fall 
within their mutual Poisson uncertainties, and there is overlap between the la confidence 
limits in their best-fit Schechter parameters (Figure IT3|) . 

We conclude that we find indications for evolution in the luminosity function over the 
2^-2^4 range, with a decrease towards lower redshifts, but only at a low significance 
level. Larger samples, such as the ones HETDEX will produce in its few firsts months of 
operation, will be required to confirm this. If real, this evolution implies that the large drop 
in the Lya luminosity function, evident at z ~ 0.3, starts to occur at z > 2. Another way 
of searching for evolution in the Lya luminosity function is to integrate it, and compare the 
implied Lya luminosity density at each redshift. This is the subject of the next section. 



7. Evolution of the Lya Luminosity Density and the Global Escape Fraction 

of Lya Photons. 

In §5.2 we measured the median escape fraction of Lya photons at 2 < z < 4 in LAEs 
to be ~ 30%. This does not represent the median escape fraction of the whole galaxy 
population at those redshifts, since LAEs will, by definition, be biased towards having high 
f esc (Lya). On the other hand, we can integrate the Lya luminosity function measured 
in the previous section to estimate the Lya luminosity density (pLya) at these redshifts. 
Comparing this observed luminosity density with that predicted from the global SFR density 
(psfr) for the entire galaxy population provides an estimate of the global escape fraction 
of Lya photons and its evolution with redshift. The above approach is equivalent to taking 
the ratio between the SFR density implied by the observed Lya luminosity density using 



Equation 2 (psFR,L ya ), and the total Psfr- This method has been applied by ICassata et al. 
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( 120 111 ). In this work we extend their a nalysis which in c luded the ICassata et al.l (1201 if ) data 
at 2 < z < 6.6, t he measurement of iGronwall et al.l ( 120071 ) at z — 3.1, and the data of 
Ouchi et al.1 (12008^ at z = 3.1, 3.7, and 5.7. We add the HETDEX Pilot Survey data 
points at 1.9 < z < 3.8, as well as the z ~ 0.3 LAE data fro m iDeharveng et al.l ( 120081 ) and 
Cowie et al.l (12010!) . the z = 2.2 data of lHayes et al.l ( 2010al). the z = 4 .5 me asurement by 
Dawson et al.l (120071). the measur ement at z — 5.7 of IShimasaku et al.l ( 120061 ). the z 
data from Kashikawa et al. ( 2006), t he data of Ouchi et al. ( 201oh at~z=6.6. and the z 



measurement of Hibon et al. 



6.5 
7.7 



A sim ilar dataset has been analyzed in this way in 



a recent submission by lHayes et al.l ( 12 01 Obi ) , although using a different set of Ha and U V 
luminosity functions at different redshifts to estimate the total SFR density. 

The top panel in Figure ITU shows psFR,L ya derived from the observed Lya luminosity 
density using Equation 2. We present our results for the full sample and for the low-z 
and high- z bins of the HETDEX Pilot Survey (red, blue, and green filled circles), as well 
as the compilation of data points calculated from the Lya luminosity functions at 0.3 < 
z < 7.7 mentioned above (black filled circles). Vertical error-bars are estimated from the 
published uncertainties in L* and 0*, and horizontal error-bars show the redshift range of the 
different samples (omitted for narrow-band surveys). A lso presented is the lat est estimate 
of the total SFR density history of the universe from iBouwens et al.l ( )2010bl ). which has 
been derived from the best to date collection of dust extinction corrected UV luminosity 
func tions at a series of red shifts between and 8, and shows a typical uncertainty of 0.17 
dex ( IBouwens et al.l 1201 Obi and reference therein). 



A source of systematic error in measuring the Lya luminosity density comes from the 
choice of the luminosity limit down to which the integration o f the luminosity fun ction is 
performed. An excellent discussion on this issue can be found in lHayes et al.l ( 12010bl ). With 
the goal of estimating the volumetric Lya escape fraction by comparison to the UV derived 
SF R density, we should i deally choose an integration limit consistent with the one used 
by IBouwens et al.l (j2010bl ) to integrate their UV luminosity functions. In this way we can 
ensure both measurements are roughly tracing the same galaxies. In the case of Lya and 
UV luminosity functions this is nontrivial, as the exact number will depend on the, mostly 
unconstrained, shape of the Lya escape fraction distribution for galaxies as a functio n UV 
luminosity. In lack of a better choice, we follow the approach of lHayes et al.l (j2010bl ). and 
integrate the Lya luminosity functions down to th e same fraction o f L* as the UV luminosity 
func tion were integrated (0.06L* =3 in the case of IBouwens et al.l (]2010bl ) ) . For consistency 
with lHayes et al.l (l2010bl ) , and in order to allow for a better comparison with their results, we 
define this limit using the lGronwall et al.l ( 120071 ) luminosity function at z = 3.1, for which the 
integration limit becomes 0.06Lq 07 = 2.66 x 10 41 erg s _1 . For all the data points in Figure 
we also shows the same measurements obtained by integrating the luminosity functions down 
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to zero luminosity (upper open circles). The unlimited integration typically overestimates 
the luminosity (SFR) density by ~ 60%. This provides a notion of the maximum impact 
that the choice of this integration limit has on the measured value of the luminosity density. 

A second source of systematic error in the above measurement comes from the role that 
IGM scattering has at reducing the observed Lya flux of sources at very high redshifts. 
Although all our previous analysis neglected the effects of IGM scattering on the Lya line, 
this approach was only well justified at our redshifts of interest (z < 4), where IGM scattering 
is negligible for our purposes (see discussion in §1 and §5). To study the escape of Lya 
photons from galaxies across a larger redshift range, we should try to incorporate the effect 
of the IGM, which is not negligible for the measurements at very high redshift (z ~ 6). As 
discussed in §1, the effects of the ISM and IGM kinematics in and aro und galaxies makes 



this correction very difficult (IDiikstra et a 



2007 



have applied a correction using the iMadaul ( 119951 ) average Lya forest opacity, and assuming 



Verhamme et al.ll2008l ). To first order, we 



that only half of the Lya line flux suffers this attenuation. The filled symbols in Figure [T 
include this correction. Raw measurements done without applying this correction are also 
shown in Figure HH as the open circles below each data point. While this correction can 
become large (~ 50%) at the highest redshifts, it impact is still within the la uncertainties 
coming from the luminosity function measurements. 

In accordance with the low significance hint of evolution presented in §6.3, in the 2 < 
z < 4 range, all the estimates of PsFR,Lya agree with each other within la. However, by 
examining the overall trend of the data points, and keeping in mind the ones at higher and 
lower redshifts, there are clearly indications for evolution in the SFR density derived from 
Lya from z ~ 7 down to z ~ 0.3, with a steady decrease towards lower redshifts across the 
2 < z < 4 range. Although the uncertainties in the 2 < z < 4 range are large, allowing any 
two datapoints to be consistent with each other, the overall trend implies a decrease in pi ya 
of roughly a factor of ~ 2 from z = 4 to 2. We stress the need for larger samples of LAEs at 
these redshifts to better constrain this evolution. 

The bottom panel of Figure [T41 shows the global average escape fraction of Lya photons, 
which is given by the ratio between psFR,Lya and Psfr at any given redshift. The average 
escape fraction derived from our Lya luminosity function over the whole 1.9 < z < 3.8 
range is (3.0li'2)%. Errors include la uncertainties in the luminosity fun c tion p arameters 
and the 0.17 dex uncertainty in the total SFR density from lBouwens et al.l ( 12010bl ). For our 
1.9 < z < 2.8 and 2.8 < z < 3.8 bins we derived a mean Lya escape fraction for the overall 
galaxy population of (2.0^9)% and (4.3+2° 2 3 )% respectively. This amount of evolution is 
not statistically significant, but we believe it to be real in the context of the overall trend 
seen in Figure HH It also does not contradict the lack of evolution in the escape fraction 
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for LAEs observed in §5.2, since, as mentioned above, the LAE selection tends to identify 
galaxies with high f esc (Lya) at any redshift, independent of the value of the escape fraction 
of the total galaxy population. The median dust extinction of a factor of ~ 3 measured in 
§4.2 implies that LAEs contribute roughly 10% of the total star formation at 2 < z < 4. 
This contribution rises to 80% by z ~ 6, implying that galaxies at these redshifts must have 
very low amounts of dust in their ISM, which is co nsistent with the very blue UV slopes of 



conti nuum selected galaxies at these high redshifts (IBouwens et al. 



2010a; 



Finke 



stein et al. 



20100 . The observed behavior is also consistent with the results of IStark et al.l (J2010|), who 
find the fraction of LBGs showing high EW Lyq em i ssion to roughly double from z = 4 to 
6. A similar result was also reported by lOuchi et al.l ( 120081 ). who measure a significant level 
of evolution in the UV luminosity function between their z = 3.7 and z = 6.6 samples of 
LAEs, which was not traced by the Lya luminosity function. 



The above escape fractions are in agreement with the result of lHayes et al.l ( 1201 Oaf ), who 
measured an overall Lya escape fraction of (5.3 ± 3.8)% at z = 2.2 by comparing the Lya 
and Ha luminosity function of narrow-band selected emission line galaxies ove r the same co 



movin g volume. On the other hand, by applying the same method used here ICassata et al. 
(120111 ) measured an average escape fraction of ~20% at z = 2.5. The difference is easily 



explained by the fact that the latter authors compared their observed Lya derived SFR 
density (which agrees with our value) to the total SFR density uncorrected by dust, which 
underestimates the true value at these redshifts. 

It is evident that a strong decrease in the Lya escape fraction of galaxies occurred 
between z ~ 6, and z ~ 2. In order to quantify this decrease we fit the data points in the 
lower panel of Figure HU using two different functional forms. First, we fit a power-law of 
the form 



log(/ MC (*)) = log(/ MC (0)) + £ log(l + z) 



(7) 



This is the same parametrization used by lHayes et al.l (l2010bl ) to fit the history of the 
global Lya escape fraction. Best fit parameters are presented in Table In order to provide 
a quantitative assessment of the impact of systematic errors in the measurement, we not only 
fit our best estimates of the escape fraction at each redshift, but also the values calculated 
ignoring the luminosity function integration limit, and the IGM correction. The best fitting 
power-laws for these three se ts of data points are shown as dotted lines in Figure O For 
comparison with the results of lHayes et al.l (j2010bf ). we should consider our raw measurement 
without including the effects of the IGM, as a correction of this type was not done in their 
work. They find best fit values of log(/ esc (0)) = —2.8 ± 0.1, and £ = 2.6 ± 0.2, in excellent 
agreement with our result. 
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While the power-law model provides a reasonable fit to the data, it seems to system- 
atically overestimate the measured values of f esc (Lya) in the 2 < z < 5 range, and under- 
estimate them in the 5 < z < 8 range. The data points in Figure O seem to indicate a 
sudden drop, or transition in the escape fraction between z = 6 and 2. A similar transition, 
in a coincident redshift range, is observed in the dust extinction derived from the UV slope 
of continuum selected galaxies (IBouwens et aJj 120091 ). Given the important role that dust 
has at regulating the escape of Lya photons, it would not be surprising if the dust content 
and the Lya escape fraction of galaxies present a similar evolution with redshift. In order 
to quantify this behavior we also fit a transition model of the following form, 



log(/ MC (z)) = ^ V " (1 - tanh(0(* - z tr )) (8) 

where z tr is the transition redshift at which the decrease in the escape fraction takes place 
{fesc = 1 for z ^> z^), and the parameter 9 determines the sharpness of the transition. Best- 
fit parameters to the measured escape fraction at each redshift, and the values without IGM 
correction, and without a luminosity function integration limit are presented in Table [3) Our 
best-fit transition model, implies a very high Lya escape fraction of ~ 80% at z ~ 6, which 
drops softly from z ~ 6 to z ~ 2, with a characteristic transition redshift at z tr = 4.0 ± 0.5, 
in order to reach a value of ~ 1% in the local universe. 

By analyzing the values in Table 3, It can be seen that, given the current uncertainties, 
the IGM correction and the choice of the luminosity function integration limit do not induce 
major changes in the best-fit parameters, especially in the case of the power-law exponent. 
The largest effect is that of the integration limit on the escape fraction at z = 0. The reason 
for this is that low L* values are measured for the Lya luminosity functions at low redshift. 
Therefore, the integration limit lays closer to L* at these redshifts, making the value of the 
luminosity density more dependent on it. 

Equation 8 predicts the average Lya luminosity of star-forming galaxies at any redshift 
given their average SFR, and it might prove useful for semi-analytic models of galaxy for- 
mation attempting to reproduce the Lya luminosity function. However, the escape fraction 
shows a very large scatter for single objects, and it might be systematically different for 
galaxy populations selected using different methods. Therefore, this relation should be used 
with caution, and only in an statistical manner. Also, this equation is only valid over the 
redshift range in which observations are available, and to the extent that current uncer- 
tainties allow. For example, given the current uncertainties, we do not consider the escape 
fraction to be properly constrained at z > 6.6. While it is tempting to interpret the slight 
drop seen in the last data point at z = 7.7 as a possible reduction in Lyct transmission due 
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to the neutralization of the IGM as we walk into the end of re-ionization, the error-bars are 
too large to allow for any significant conclusions. 

8. Conclusions 

For a sample of LAEs at 1.9 < z < 3.8, detected by means of blind integral field spec- 
troscopy of blank extragalactic fields having deep broad-band optical imaging, we were able 
to measure the basic quantities SFR, E(B — V), UV luminosity, Lya EW, and f esc (Lya). 
From these quantities and the correlations observed between them we conclude: 

• Over the 2 < z < 4 range LAEs show no evolution in the average dust content of 
their ISM, parameterized by the dust reddening E(B — V), and measured from the UV 
continuum slope. These objects show a mean {E(B — V)) = 0.13 ±0.01, implying that 
dust absorbs ~ 70% of the UV photons produced in these galaxies. While one third 
of LAEs down to our luminosity limit are consistent with being dust-free, the level of 
dust extinction measured for the rest of the sample is significant, and should not be 
neglected. 

• At z < 3, we see the possible appearance of a rare (6/89 objects) population of highly 
reddened {E(B — V) > 0.45) LAEs, typically showing high Lya EWs. Two of these 
objects are in the MUNICS field where we do not have deep X-ray data to reject 
AGNs from our sample. The remaining four objects show low association probabilities 
(< 0.6) with their broad-band counterparts, casting doubt on the validity of our UV 
slope and EW measurements. The presence of these objects in the sample does not 
affect significantly the average dust properties of LAEs at the low redshift end of our 
range. If real, these objects are of great interest since their presence could indicate that 
the fraction of dusty LAEs grows towards lower redshifts. Followup of these objects is 
necessary to confirm this. 

• The Lya EWs of LAEs in our sample are consistent with the expectations for nor- 
mal stellar populations with metallicities within 1/50 Z Q and solar. We do not find 
it necessary to invoke a top-heavy IMF, the presence of population III stars, or an 
enhancement of the EW due to a clumpy dust distribution in a multi-phase ISM, to 
explain the observed EWs. 

• LAEs in our sample show a median uncorrected UV derived SFR ~ 11 M Q yr -1 . 
Correcting the UV luminosities for dust extinction increases this median value to 
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SFR ~ 34 M yr _1 , implying that assuming LAEs to be dust-free galaxies can trans- 
late into large underestimates of their SFRs. The ratio between the observed (i.e. 
uncorrected for dust) UV and Lya derived SFRs shows a median value of 0.83. Nei- 
ther this ratio, nor the Lya EW, show significant evolution with redshift across the 
2 < z < 4 range. These two quantities also show no dependence with E(B — V), 
implying that the ratio between the amount of dust extinction seen by Lya photons 
and that seen by UV photons is independent of the dust-content of the galaxies' ISM. 
This finding is at odds with the expectation of models in which a clumpy distribution 
of dust in a multi-phase ISM promotes the escape of Lya photons over that of UV 
continuum photons. It also implies that some combination of ISM geometry and kine- 
matics reduces the amount of extinction seen by Lya photons from that expected in 
a static and homogeneous ISM, but only up to the point where it is similar to that 
experienced by continuum photons. 

• The escape fraction of Lya photons from LAEs, given by the ratio between the observed 
Lya luminosity and that predicted from the dust-corrected UV SFR, shows a median 
value of 29%. A large scatter is seen around this number, with objects ranging from 
a few percent to 100%. Both the median value, and the range of observed escape 
fractions in LAEs, show no evolution across the 2 < z < 4 range sampled by our 
objects, and does not seem to evolve all the way down to z — 0.3. Since LAE selection 
is biased to include objects of high escape fractions at any combination of dust content, 
redshift and survey limiting luminosity, it is not surprising that this parameter shows 
little or no evolution. This lack of evolution in LAEs does not imply that the Lya 
escape fraction for the overall galaxy population is not evolving. 

• The Lya escape fraction of LAEs shows a clear correlation with E(B — V). The 
correlation follows what is expected for a value of q = 1, where q is the ratio between 
the dust opacity seen by Lya and that seen by continuum photons. This behavior is 
consistent with what is observed for LAEs at z — 0.3, implying that not only the value 
of the escape fraction, but also its dependence with dust content, do not evolve with 
redshift. While other galaxies not identified by the LAE selection can have q > 1, and 
show low Lya EWs, lack of Lya emission, and even Lya in absorption, the lack of 
objects with q <C 1 confirms that preferential escape of Lya photons over continuum 
photons in the presence of a clumpy dust distribution is not a common process in 
galaxies. 

We also measure the Lya luminosity function across our redshift range. Integration of 
the luminosity function yields a measurement of the Lya luminosity density in our sampled 
volume. Using our data, and a compilation of measurements of the Lya luminosity function 
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at different redshifts from the literature, we are able to trace the evolution of p Lya with 
redshift from z = 0.3 to z = 7.7. Comparing the observed value of pL ya with the expected 
density implied by the SFR history of the universe, allows a measurement of the evolution 
of the average Lya escape fraction for the overall galaxy population in this redshift range. 
From these measurements we conclude the following: 



The observed Lya luminosity function is well matched to previous measurements in 
the literature, especially in the L(Lya) < 10 43 erg s -1 range typically measured by 
previous studies. Given the large volume sampled by the HETDEX Pilot Survey, we 
are able to find many objects above this luminosity. Both the redshift distribution 
and the luminosity function show hints of a decrease in the number density of LAEs 
of roughly a factor of 2 from z = 4 to 2, although this decrease is not statistically 
significant and larger samples are required before it can be confirmed. In any case, 
this decrease goes in the right direction and is consistent with what is expected from 
the observed drop in the Lya escape fraction for the overall galaxy population. 

The Lya luminosity density is observed to increase steadily from z = 0.3to;?~2, 
following the behavior of the SFR history of the universe. Over this range, the average 
Lya escape fraction increases very slowly from ~ 1% to ~ 5%. At z > 2 the increase 
in p Lya starts to flatten, and a decline is observed around z ~ 6. This behavior 
is accompanied by a decrease in psfr immediately after z = 2, implying that over 
the 2 < z < 6 range, the average Lya escape fraction in galaxies increases steadily 
from the ~ 5% up to ~ 80% by z = 6. Current measurements of the luminosity 
function at higher redshifts do not constrain the behavior of the escape fraction beyond 
z = 6.6. This drop in the average escape fraction of Lya photons with cosmic time 
is consistent with the increase in the dust-content of star forming galaxies, which is 
expected from the chemical enrichment of these objects as star formation proceeds, 
and is ob served as a reddening i n the UV slope of star forming galaxies towards lower 



redshifts (IBouwens et al.ll2010al ; iFinkelstein et al.ll2010h 



Equation 8 provides a useful analytical form which describes the history of the average 
Lya escape fraction of galaxies as a function of redshift. This equation can prove useful 
to predict the expected average Lya luminosity of galaxies in numerical simulations 
and semi-analytical models. The reader must keep in mind that galaxies do not show 
a single value of the escape fraction at any given redshift, but rather a relatively 
broad (and mostly unconstrained) distribution, so this equation can only be used in a 
statistical sense. It must also be kept in mind that the behavior of the escape fraction 
at z > 6.6 is still unconstrained. 
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These last few points have interesting consequences regarding the potential that observa- 
tions of LAEs at very high redshifts (z > 7) have to detect the effects of cosmic re-ionization. 
Our results imply that at these redshifts, dust is no longer an important factor setting the 
average escape fraction of Lyct photons in galaxies. Therefore, a significant drop in the Lya 
escape fraction could be more easily interpreted as being caused by an increased neutral 
fraction in the IGM. 
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Fig. 1. — Limiting Lya luminosity (5a) as a function of redshift for the survey. The survey 
depth varies across the observed area due to changes in atmospheric transparency, Galactic 
extinction, and instrumental configuration. Hence, the background color indicates the frac- 
tion of the total survey area over which a given limit is reached. White points mark the 
redshift and Lya luminosities (with error-bars) of the 99 objects classified as LAEs. The 
dotted black and white lines show the mean and best limits over the whole survey respec- 
tively. Even below this last limit, the completeness of the survey is not zero, explaining why 
we see 2 points below this curve. 
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Fig. 2. — Redshift distribution of the 99 LAEs in the Pilot Survey (solid histogram). 
Error-bars represent Poisson uncertainties only. Also shown is the incompleteness-corrected 
predicted redshift di stribution (dotted line) given by our flux limit and assuming the 
Gronwall et all (120071 ) Lya luminosity function with no evolution over 2 < z < 4. 
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Fig. 3. — UV continuum slope as a function of redshift for the 89 LAEs with broad-band 
optical counterparts . Objects are color co ded by field. The right axis shows the equivalent 
E(B-V) assuming a ICalzetti et al.l (120001 ) attenuation law. The horizontal lines mark the 
assumed intrinsic UV slope corresponding to a dust-free stellar population (/3 = —2.23, 
solid line), and the mean for the whole sample (dotted line). Also shown are the mean UV 
slopes for two redshift bins at 1.9 < z < 2.8 and 2.8 < z < 3.8 (black squares), with two sets 
of error-bars corresponding to the standard deviation in /3 within each bin (large error-bars) 
and the formal error in the mean (small error-bars). 
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Fig. 4. — E(B — V) distributi on of LAE s in ou r sa mple (Poisson error- bars), together with 
that of BX/LBGs taken from ferb et all kood \ and iReddv et all J2008h (solid histograms). 
The median of each distribution is marked by the vertical dashed lines. 
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Fig. 5. — Rest-frame Lya EW distribution of LAEs in our sample (dashed black histogram). 
The distributions for low (E(B — V) < 0.45) and high (E(B — V) > 0.45) reddening objects 
are shown (blue and red histograms respectively). Also shown are the best-fit exponential 
distribution (N oc exp [—EW/wq]) to the whole sample (wq = 77 ± 7A, solid black line) and 
the low reddening sample (wq = 74 ± 7 A, dotted blue line). 
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Fig. 6. — UV versus Lya derived SFRs for the 83 LAEs in the final sample. Values are 
not corrected for dust extinction. The solid line shows the median SFR(Lya) to SFR(UV) 
ratio of 0.83. The expected range for dust-free normal stellar populations is marked by the 
dashed lines. Dotted lines mark ratios of 0.01, 0.1, 1, 10, and 100. 
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Fig. 7. — Rest-frame Lya and SFR(Lya) to SFR{UV) ratio (not corrected for dust) 
as a function of redshift. The median EW of 71A and ratio of 0.83 are marked by solid 
horizontal lines. The dotted lines on the top panel indicate the maximum EW range for 
young normal stellar populations with metallicities between solar and one 1/50 solar from 
Schaererl (120031 ). Dotted lines in the bottom panel display the allowed range in the SFR(Lya) 
to SFR(UV) ratio for dust-free normal stellar populations. The open boxes show the median 
EW and ratio for the two redshift bins at 1.9 < z < 2.8 and 2.8 < z < 3.8. 
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Fig. 8. — Rest-frame Lya EW and SFR(Lya) to SFR(UV) ratio (not corrected for dust) 
as a function of E(B-V). Symbols are the same as in Figure [3 
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Fig. 9. — Same as Figure El but with SFR(UV) corrected for dust. Error-bars include the 
uncertainty in the correction. The solid line marks the median escape fraction of 29%. 
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Fig. 10. — Escape fraction of Lyct photons as a function of redshift for the 83 LAEs in the 
final sample. The solid horizontal line denotes the median escape fraction of 29%. Also shown 
is the median escape fraction for the two redshift bins at 1.9 < z < 2.8 and 2.8 < z < 3.8 
(open red stars), with error-bars corresponding to the standard deviation of logf/^) within 
each bin. The escape fractions of LAEs at z = 0.3 with their median from lAtek et al.l ( 2009 ) 
(green triangles, red open square) are also displayed. 
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Fig. 11. — Lya escape fraction as a function of E(B-V). Dashed lines show the expected 
correlation for different values of the parameter q = TT, vn ./Tx=^w- Th e red line displays the 
relation for LBGs showing Lya in emmission fromlKornei et al.l (120101 ) . Green triangles show 
the values for z ~ 0.3 LAEs from ktek et al.l (|2Q09h . 
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Fig. 12. — Lya luminosity function of the HETDEX Pilot Survey sample of 80 LAEs 
in COSMOS and HDF-N, shown before and after applying the completeness correction 
(open black and filled red circles respectively). Poisson error-bars are included. Also 
displayed are the completeness corrected luminosity function for the two redshift bins at 
1.9 < z < 2.8 and 2.8 < z < 3.8 (bl u e and green star s respectively), and t he luminosity func 



tions of Ivan Breukelen et all (120051) ; iGronwall et all ( 120071 ) ; lOuchi et all ( 120081 ) ; lHayes et al. 



( 12010a) ) .and ICassata et all ( 120111 ). Schechter fits to the full sample, as well as the low-z and 
high-z samples, are also presented (solid red, blue, and green curves respectively). The red 
dashed line denotes the best Schechter fit to the L(Lya) < 10 43 erg s _1 bins. 
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Fig. 13. — Contours show 1, 2, and 3a confidence limits for the luminosity function param- 
eters L* and (f>* . Stars show our results for the full sample and the two redshift bins a t 
1.9 < z < 2.8 a n d 2.8 < z < 3.8. The parameters estima t ed by van Breukelen et al . (2005); 



Gronwall et all (120071 ): lOuchi et all (120081 ): lHaves et all (I2010af ). and 
are also presented (filled circles). 
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Fig. 14. — Top panel: SFR density (psfr) as a function of redshift. The solid and dot- 



ted lines show the total Psfr from iBouwens et al.l (l2010bl ) and its typical uncertainty of 
0.17 dex. Blue, green, and red filled circles show PsFR,Lya derived from the Lya luminosity 
function in the two redshift bins at 1.9 < z < 2.8 and 2.8 < z < 3.8, as well as for the 
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shown by the open circles below each measurement. Values computed integrating the Lya 
luminosity functions all the way down to L(Lya) = are shown by the open circles above 
each measurement. Bottom panel: Escape fraction of Lya photons for the overall galaxy 
population, derived from the ratio between the Lya derived psFR,L ya and the total value at 
each redshift. The dashed line marks an escape fraction of 100%. Solid lines shows our best 
fit to the data given by Equation 6, while dotted lines show the best fit powerlaw functions. 
Purple, orange, and cyan colors indicate fits to the escape fraction measurements including 
an IGM correction and an integration limit for the luminosity function, ignoring the IGM 
correction, and ignoring the l uminosity function in tegration limit respectively. The black 
dashed line shows the result of lHayes et al.l ( 12010bl ). 
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Table I. Properties of HETDEX Pilot Survey LAEs 



ID a 



L(Lya) 



L » b 

1/.1500A 



E(B-V) fesc(Lya) EW (Lya) 









10 42 erg s _1 


10 28 erg s _1 Hz -1 






map" 






A 


HPS- 3 


3 


09 


14.4±2.8 


12.1±1.6 


-0.9±0.4 





,27±0 


.08 


06+2?, 8 , 


734-16 


HPS-6 


2. 


78 


20.1±2.2 


19.5±1.6 


-1.4±0.2 


o 


18±0 


.06 


12+°/i° 
u,± — 0.06 


58±8 


HPS- 11 


2 


78 


11.5±2.5 


18.7±1.1 


-2.2±0.2 


o 


00±0 


.05 


40+S - 2 1 


28±6 


HPS- 13 


3 


.32 


10.1±2.0 


9.3±1.9 


-1.2±0.5 


o 


,20±0 


.11 


10+SiI 

u -- LU -0.07 


62±16 


HPS-17 


2 


78 


6.9±1.9 


13.2±1.5 


-2.5±0.4 


o 


00±0 


.08 


34+° f 3 


23±6 


HPS-22 


2. 


.77 


9.8±2.7 


1.9±1.0 


-0.6±1.2 


o 


33±0 


.25 


14+i?5 

u- J-t_0.14 


3404-187 


HPS-25 


2 


.55 


32.1±4.7 


9.1±1.2 


-0.2±0.3 





,41±0 


.08 


05 +0 ,'°,5 

u-uu — 0.03 


252±50 


HPS- 34 


2, 


.76 


11.2±2.8 


11.7±1.0 


-1.9±0.2 


o 


,07±0 


.06 


u "°° — 0.16 


47±12 


HPS-51 


3 


10 


5.9±2.9 


22.4±2.6 


-1.0±0.3 


o 


26±0 


.08 


u-u — 0.01 


16±8 


HPS-53 


3 


.57 


13.0±4.2 
















HPS-62 


2, 


08 


17.1±5.5 


6.9±0.8 


-1.3±0.3 





,18±0 


.07 


0.29 + °?2 

" — 0.17 


139±49 


HPS-82 


2, 


25 


29.1±7.7 


1.4±0.7 


1.7±1.4 





,79±0 


.29 


0.01+°" 

— 0.01 


2213±1548 


HPS- 84 


3 


25 


24.3±6.1 


15.1±3.8 


-1.3±0.5 





19±0 


.11 


0.17+°'?? 

— 0.12 


91±30 


HPS- 89 


2, 


54 


14.4±3.0 
















HPS-91 


3 


.00 


10.3+.3.4 


23.7±3.2 


-0.8±0.3 





.28±0 


.08 


0.02+0- 02 

— 0.01 


27±9 


HPS-92 


3 


.67 


13.6±4.4 


44.9±5.7 


-1.6±0.3 





,14±0 


.07 


0.05+°°6 

— 0.O3 


16±5 


HPS-93 


2 


26 


20.9+.5.0 














HPS-95 


2 


.45 


13.4±4.8 


2.1±0.6 


-1.8±1.2 





09±0 


.25 


1.72+J™ 9 


322±155 


HPS-99 


3. 


.01 


22.8±4.7 


6.0±2.2 


-0.4±0.7 





37±0. 


.16 


0.07+°^ 

— O.Oo 


258+.102 


HPS- 109 


3 


.21 


22.5±5.9 


47.5±5.6 


-0.9±0.2 





28±0 


.06 


0.02+°-°? 

— 0.01 


29±8 


HPS-111 


3 


.18 


11.2±3.7 


24.0+.2.9 


-1.7±0.3 





,11±0 


.07 


0-10+n'ofi 

— ().(.)() 


24±8 


HPS-124 


3 


.74 


13.0+.6.3 


10.1+.3.7 


-3.0±0.9 





00±0 


.19 


0.83+"? 

— . 5 1 


51±27 


HPS-126 


2. 


.83 


106.7+.9.1 


3.5±3.2 


1.9±1.4 





,82±0 


.28 


0.01+"'"? 

— 0.01 


3338±3038 


HPS-127 


2. 


.54 


9.0±3.6 


10.2±0.8 


-1.5±0.3 





15±0. 


.08 


0.14+"-": 


48±19 


HPS- 142 


2, 


58 


9.1±2.2 


20.2±1.1 


-1.1±0.2 





22±0. 


.05 


0.04+°-°? 

— 0.1)2 


26±6 


HPS- 144 


2 


73 


2.7±1.3 


1.0±0.5 


1.3±1.0 





70±0 


.20 


O.OO+Ho 

— 0.00 


270±187 


HPS- 145 


2 


.18 


26.5±3.6 


5.4±0.5 


0.1±0.3 





.48±0 


.07 


0.03+°°? 

— 0.02 


380±66 


HPS-150 


2 


.90 


18.1±4.2 


17.7+.1.2 


-1.5±0.2 





15±0 


.05 


0.16t°o7 


55±13 


HPS-153 


2 


.71 


16.3+.3.2 


5.1±1.0 


-0.9±0.4 





26±0. 


.08 


n 18+022 


198±50 


HPS-154 


2 


.87 


6.2±1.6 


2.5±1.0 


-1.1±0.9 





,22±0 


.19 


19+ 103 


148±64 


HPS- 160 


2. 


.43 


6.6±3.2 


0.4±0.4 


0.1±2.2 





46±0. 


.44 


14+ 903 


1306±1732 


HPS-161 


3 


25 


35.1+.3.7 


31.6+.2.4 


-0.4±0.2 





,37±0 


.05 


02+ 01 
u - uz -o.oi 


76±9 


HPS- 164 


2. 


.45 


10.1+.5.3 


7.4±0.8 


-1.3±0.2 





20±0. 


.06 




79±42 


HPS- 168 


3 


.45 


36.4±3.0 


7.5±1.1 


-2.0±0.3 





04±0 


.07 


2.08tg-j° 


238±35 


HPS-174 


3. 


.45 


2.7±2.0 


2.1±0.9 


-2.5±1.1 





00±0 


.22 




58±45 


HPS- 182 


2. 


.43 


10.4±2.2 


4.5±0.4 


-2.0±0.3 





04±0. 


.08 


99+ 1 11 


114±27 


HPS-183 


2 


.16 


8.6±5.3 


3.2±0.4 


-2.1±0.4 





03±0. 


.09 


1 - OD -0.90 


128±82 


HPS- 184 


3 


.21 


4.4±3.0 


4.4±1.5 


-1.8±0.9 





08±0 


.19 


29 +1 ' 53 
u '^ a -0.26 


51±36 


HPS- 189 


2. 


.45 


4.9±2.9 


4.5±0.6 


-2.0±0.3 





05±0 


.07 


43+0 50 


54±32 


HPS- 190 


2 


.28 


6.0±1.6 
















HPS- 194 


2 


.29 


23.5+.1.8 


10.6+.0.8 


-1.8±0.2 





09±0 


.06 


62+ 44 
U - DZ ~0.26 


114±13 


HPS-196 


2. 


.65 


12.3+.2.0 


7.4±1.0 


0.4±0.3 





53±0. 


.07 


o m +0 01 
u - ul -o.oo 


134±27 


HPS-197 


2 


.44 


7.1±2.6 


2.2±1.0 


-1.9±1.8 





,07±0 


.36 


±.U9_ .7 


160±93 


HPS-205 


2 


91 


12.7+.3.5 


2.1±0.9 


-1.0±0.9 





26±0 


.19 


34+ 1 - 71 
u " 34 -0.30 


372±174 
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Table 1 — Continued 



ID a 


z 


L(Lya) 


iaISOOA 




P 




E(B - V) 


fesc(Lya) 


EW (Lya 


HPS- 


■207 


2. 


.71 


5.0±1.7 


2.1±0.8 


-2 


,9±1 


.2 


0. 


00±0. 


.25 




97±46 


HPS- 


■210 


3 


.49 


9.5±3.0 


9.6±2.4 


-1 


,3±0 


.5 





18±0, 


.11 


1 I +0.19 
u - il -0.08 


56±20 


HPS- 


■213 


3 


.30 


11.0±2.8 


11.6±1.3 


-0 


.5±0 


3 





35±0. 


.07 


02+ 02 
u ' uz -0.01 


63+17 


HPS- 


-214 


3 


.30 


6.6±3.1 


1.4±0.5 


-2 


,7±1 


.2 





00±0 


24 


o 1 1+28.94 


202+102 


HPS- 


■223 


2. 


.31 


12.9+3.5 


2.0±0.5 


-1 


,2±0 


.6 


0. 


21±0. 


13 


55+ 1 ' 41 
u.oj_q 42 


373+142 


HPS- 


-229 


3 


.04 


31.6±3.5 


30.5±1.9 


-1 


,6±0 


.2 





14±0 


05 


0.l8t°;J 2 7 


55+6 


HPS- 


■231 


2 


.72 


16.1±4.1 


1.8±0.5 


-1 


,9±0 


.8 





07±0. 


.18 


2.95t 1 1 2 7 P 


459+190 


HPS- 


■244 


2. 


.10 


2.6±1.2 


1.7±0.4 


-2 


.3±0 


.7 





00±0 


.15 


i nn+3- 04 


71+38 


HPS- 


■249 


3 


.27 


5.7+2.2 


2.5±0.8 


-2 


,6±0 


.7 


0. 


00±0 


.14 


lA8+_tf 5 


98+44 


HPS- 


■251 


2 


.29 


14.3±4.0 


5.2±0.5 


-1 


.9±0 


.3 





07±0. 


.08 


0.88t°; 9 5 f 


140+43 


HPS- 


■253 


3 


.18 


15.4+3.0 


12.9±1.4 


-1 


,7±0. 


.2 





10±0. 


.06 


n 2q +0.24 


62+13 


HPS- 


-256 


2 


.49 


13.9+3.5 


3.5±0.6 


-1 


,7±0 


.4 





10±0. 


.09 


1 01 + 1 ' 49 


206+65 


HPS- 


■258 


2 


.81 


19.3+2.4 


13.8±0.9 


-0 


.8±0. 


.2 


0. 


28±0 


.06 


06+ 05 
u - uo -0.03 


88+12 


HPS- 


-263 


2 


.43 


9.2+3.0 


9.7±0.6 


-1 


.7±0 


.2 





10±0. 


06 


23 +0 18 


49+16 


HPS- 


■266 


2 


.20 


13.8+1.9 


- 




- 










- 


- 


HPS- 


-269 


2 


.57 


6.2+1.6 


3.6±0.5 


-1 


.9±0 


.5 





07±0. 


.10 


n .-+0.98 


87+26 


HPS- 


■273 


3 


.64 


15.0+5.7 


- 




- 










- 


- 


HPS- 


■274 


2 


.87 


10.7+2.0 


9.8±0.9 


-1 


,3±0 


.2 





19±0. 


.06 


1 1 + ' 08 
U - ii -0.05 


62+12 


HPS- 


-283 


3 


.30 


19.3+2.8 


14.4±1.6 


-1 


,3±0 


.2 





20±0. 


.06 


O 14+ 011 


76+13 


HPS- 


■286 


2 


.23 


9.2+2.0 


7.6±1.5 


-2 


1±0 


.5 





03±0. 


.11 


61+ 1 ' 13 


59+18 


HPS- 


■287 


3 


.32 


4.7+2.1 


2.4±0.9 


-1 


,5±1. 


.2 





14±0. 


.24 


n oq+2.88 
u.oo_ 2g 


107+55 


HPS- 


■288 


3 


.04 


8.4±2.1 


12. Oil. 2 


-1 


,8±0. 


.3 





10±0 


.06 


n 10+O.I6 


36+9 


HPS- 


■292 


2 


.87 


19.6+2.6 


6.3±0.8 


-1 


,6±0 


.4 





12±0 


.08 


63+ ' 75 


166+30 


HPS- 


-296 


2 


.84 


5.8+2.2 


6.3±1.0 


-1 


,6±0 


.4 





13±0 


.10 


17+ - 28 


49+19 


HPS- 


■306 


2 


.44 


14.8+2.9 


11.8±0.8 


-1 


,9±0 


.2 


0. 


07±0. 


.05 


43+ - 29 


62+13 


HPS- 


-310 


3 


.07 


7.6±1.9 


7.2±1.2 


-1 


,5±0 


.4 


0. 


15±0 


09 


16+ ' 21 
u - lu -o.io 


58+16 


HPS- 


■313 


2 


.10 


6.7+3.0 


24.9±0.8 


-1 


.5±0 


1 





15±0 


.04 


04+ 03 


14+6 


HPS- 


■314 


2 


.63 


6.9+2.1 


- 




- 










- 


- 


HPS- 


■315 


3 


.07 


5.9±1.8 


14.7±1.5 


-1 


,5±0 


.3 





14±0 


07 


n n? +0.06 


21+6 


HPS- 


■316 


2. 


.81 


13.1+3.4 


13.2±1.1 


-2 


,0±0 


.2 


0. 


04±0. 


06 


n 44+ - 36 

u - 44 -0.18 


48+13 


HPS- 


-318 


2 


.46 


11.6+3.8 


13.4±0.8 


-1 


,3±0 


1 





,18±0, 


.05 


o.iot8:& 


49+16 


HPS- 


-327 


2 


.25 


4.7±1.9 




















HPS- 


-338 


2 


.60 


15.2+4.0 


1.7±0.9 


-1 


.9±1 


.1 





06±0. 


.23 


3.30lf 9 f 


452+260 


HPS- 


■341 


2 


.93 


8.4+2.3 


8.0±1.2 


-2 


,2±0 


.5 





01±0 


.10 


0.6llJ;88 


50+15 


HPS- 


■360 


2 


.92 


11.5+3.0 


7.1±2.0 


-1 


,3±0 


.5 





19±0 


12 


o.i8i8:?5 


91+33 


HPS- 


■370 


3 


.18 


8.7+2.5 


5.2±1.3 


-2 


,0±0 


.6 


0. 


04±0. 


12 


75+ 1 61 
u - '°-0.34 


81+28 


HPS- 


-372 


2 


.76 


5.5±1.4 


1.4±1.3 


-1 


.9±2 


.5 


0. 


07±0. 


.51 


-i oq+165.25 
± - ,3z -0.94 


194+183 


HPS- 


■373 


2 


.91 


11.3+2.6 




















HPS- 


■389 


2 


.59 


10.2+1.9 


7.9±1.1 


-1 


.5±0 


3 





14±0 


.08 


21 + ' 25 


70+16 


HPS- 


■391 


2. 


.96 


17.4±4.1 


8.4±2.0 


-1 


,6±0 


.5 





12±0 


.10 


44+°' 72 
u - 44 -0.29 


110+34 


HPS- 


■395 


2 


.27 


6.6+2.8 


10.2±1.2 


-1 


.9±0 


.3 


0. 


07±0. 


.07 


0-22l8;?4 


32+14 


HPS- 


■402 


2 


.97 


11.2±1.6 


5.1±1.2 


-2 


.1±0. 


.6 





02±0. 


.13 


1 on -1 " 2 ' 78 


105+28 


HPS- 


■403 


3 


.18 


7.5±1.6 


9.6±1.9 


-1 


,9±0 


.4 





08±0. 


.10 


25+ ' 38 


39+10 


HPS- 


■415 


3 


.37 


10.5+2.5 


6.1±1.2 


-2 


,0±0 


.5 





06±0. 


.11 


65+ 1 ' 19 


86+25 


HPS- 


■419 


2 


.24 


8.1±1.4 


6.0±0.8 


-1 


,7±0. 


.4 





12±0 


.08 


0.29l°;? 7 


71+15 


HPS- 


-420 


2 


.93 


12.1±2.5 


5.5±1.7 


-1 


,3±0 


.6 





19±0. 


13 


r\ or)+0.52 
u - zz -0.16 


125+47 



Table 1 — Continued 



ID* z L(Lya) ^ 1500 / P E(B-V) fesc(Lya) EW (Lya) 



HPS- 


-426 


3.41 


6.6±1.6 


6.4±1.5 


-1.3±0.5 


0.18±0.10 




58±18 


HPS- 


-428 


3.34 


13.0±2.3 


22.0±2.3 


-1.4±0.3 


0.16±0.07 


o.o8ir 8 4 


32±16 


HPS- 


-434 


2.27 


3.9±1.2 


1.0±0.4 


-2.5±1.6 


0.00±0.32 


o fjc+51 .32 
z - DO -1.36 


180±139 


HPS- 


-436 


2.42 


2.7±1.0 


4.1±0.7 


-2.8±0.6 


0.00±0.12 


42+ 94 


27±10 


HPS- 


-447 


3.13 


5.0±1.1 


14.8±2.1 


-1.6±0.3 


0.13±0.08 


o.o6±r i 


17±4 


HPS- 


-462 


2.21 


27.4±2.9 


8.3±0.8 


-1.8±0.3 


0.08±0.08 


0.98±J;J§ 


169±27 


HPS- 


-466 


3.24 


18.2±2.1 


31.3±3.4 


-1.5±0.2 


0.15±0.06 


09+ - 07 
u - uy -0.04 


32±4 


HPS 


-467 


2.80 


5.0±1.8 












HPS 


-474 


2.28 


4.3±2.4 


3.9±0.4 


-1.9±0.3 


0.07±0.07 


36+ 39 


56±32 



a ID corresponds to that in Table 3 of lAdams et alj d201lh . Equatorial coordinates and line fluxes 
are provided there. 

b Dashes indicate objects with no broad-band counterpart. 



Table 2. Lya luminosity function Best Fit Schechter Parameters, Luminosity and SFR 

Density 



Sample 


a a 


4>* 


L* 


PLya 


PSFR,Lya 






10~ 4 Mpc- 3 


10 43 erg s _1 


10 39 erg s^ 1 Mpc -3 


10 _3 M Q yr" 1 Mpc" 3 


1.9 < z < 3.8 
1.9 < z < 3.8 
1.9 < z < 3.8, L(Lya) < 10 43 
1.9 < z < 2.8 
2.8 < z < 3.8 


-1.7 

-1.5 
-1.7 
-1.7 
-1.7 


r, r, + 3.9 

z - z -1.3 

2.9T 4 ; 4 
6.7li° 9 6 

i.o±S:S 

o e+28.3 
z - D -2.2 


L - z <->_ 52 

i m +0 ' 67 

J - ul -0.41 

60+ 2 " 

1 R3+ 9 ' 46 
± - 1± -0.74 


5 1+ 2 ' 5 
°- 1 -1.6 

4.3l 2 ; 9 
6.8±I;f 

n 4 +2.7 

e c + 12.0 
°-°-2.6 


4.6t 2 ' 2 
3-9±iii 
6.2+J; 9 

5.0l 2 ° 3 9 



a Fixcd parameter 



Table 3. Lya Escape Fraction History Best Fit Paramenters 



Function 


Data Points 


log(/ esc (0)) 




Ztr 


e 


X 2 




IGM corr + LF limit 


-2.7 ±0.2 


2.4 ±0.3 






i.i 


Power Law 


No IGM corr 


-2.7 ±0.2 


2.2 ±0.3 






1.0 




No LF limit 


-2.4 ±0.2 


2.2 ±0.3 






1.2 




IGM corr + LF limit 


-2.1 ± 0.3 




4.0 ±0.5 


0.4 ±0.1 


0.41 


Transition 


No IGM corr 


-2.2 ± 0.3 




4.3 ±0.6 


0.4 ±0.1 


0.38 




No LF limit 


-1.7 ±0.2 




4.1 ±0.4 


0.5 ± 0.2 


0.39 



